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ABSTRACT 

We present three-dimensional kinematics of Sagittarius (Sgr) trailing tidal debris in six fields located 70- 
130° along the stream from the Sgr dwarf galaxy core. The data are from our proper-motion (PM) survey of 
Kapteyn's Selected Areas, in which we have measured accurate PMs to faint magnitudes in ~ 40' x 40' fields 
evenly spaced across the sky. The radial velocity (RV) signature of Sgr has been identified among our follow- 
up spectroscopic data in four of the six fields and combined with mean PMs of spectroscopically-confirmed 
members to derive space motions of Sgr debris based on ~ 15-64 confirmed stream members per field. These 
kinematics are compared to predictions of the iLaw & Majewskil (l2010a) model of Sgr disruption; we find 
reasonable agreement with model predictions in RVs and PMs along Galactic latitude. However, an upward 
adjustment of the Local Standard of Rest velocity (0lsr) from its standard 220 km s"' to at least 232 ± 14 
km s"' (and possibly as high as 264 ±23 km s"') is necessary to bring 3-D model debris kinematics and our 
measurements into agreement. Satisfactory model fits that simultaneously reproduce known position, distance, 
and radial velocity trends of the Sgr tidal streams, while significantly increasing 0lsr, could only be achieved 
by in creasing the Galactic bulg e and disk mass while leaving the dark matter halo fixed to the best-fit values 
from lLaw & Majewskil (l2010al) . We derive low-resolution spectroscopic abundances along this stretch of the 
Sgr stream and find a constant [Fe/H] ~ -1.15 (with ^ 0.5 dex scatter in each field - typical for dwarf galaxy 
populations) among the four fields with reliable measurements. A constant metallicity suggests that debris 
along the ^ 60° span of this study was all stripped from Sgr on the same orbital passage. 

Subject headings: Galaxies: individual: (Sagittarius dwarf spheroidal) — Galaxy: fundamental parameters — 
Galaxy: kinematics and dynamics — Galaxy: structure 



1. INTRODUCTION 

With the profusion of data provided in recent years by 
deep, large-area photometric surveys such as the Two Mi- 
cron All Sky Survey (2MASS) and Sloan Digital Sky Sur- 
vey (SDSS), a wealth of stellar substructure has been uncov- 
ered in the Milky Way (MW) halo. The finding and subse- 
quent mapping of numerou s stellar tidal streams and overden- 
sities (e.g., Sagittarius — Ibata et aTl 120011; iMaiewski et ajJ 
2003; Belokurov et al. 2006; Monocer os — Newberg et aTl 
20021; llbataetal. 2003; Yannv et aT ^OO l od i er SDSS 
streani s — 'Grillmair 2009; Belokurov e t al.ll2007l; iGrillmaiii 
l2006alfbl iGr illmair & Dionatos 2006) h as borne out the idea 
jSearle &Zina,1978; ,Maiewski.,1993i: iMaiewski et al.lll996l) 
that remnants of accreted dwarf galaxies make up much of 
the stellar halo of the Milky Way. The direct confirma- 
tion of the accretion of late-infalling subhalos via discov- 
ery of ubiquitous long-lived, coherent tidal debris streams 
has provided strong constraints on models of small-scale hi- 
erarchical structure formation under the prevailing A-Cold 
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Dark Matter (ACDM) co s mology (e . g.. lAbadi et a lj | 200 3|; 
iBuUock & JohnstonI l200l iFont et all [20061) . Furthermore, 
because the tidal streams retain the kinematical signatures 
of the orbits of their progenitors (i.e., angular momentum 
and energy), stellar debris in the streams can be used as 
sensitive pro bes of the unde rlying Galactic gravitational po- 
tential (e;g;J|ohnston_eta^ Ibata etal. 2001; HelnJ 
20041; iMart mez-Delgado et all ,2004; .Johnston et al., .20051 
Lawetalj |2005; MajewskietaiJlQOl). 

The best-known and only widely agreed-upon case of 
a presently visible dwarf galaxy undergoing tidal disrup- 
tion in the Milky Way halo is the Sagittarius (Sgr) dwarf 
spheroidal (dSphfl The cor e of this galaxy was first dis- 
covered by llbata et alj (119941) in a kinematical study of the 
outer Galactic bulge, with the first large-scale mapping of the 
Sgr le ading and trailing tidal arms done by Majewsk i et alJ 
(120031) using 2MASS M-giant stars. Various studies 
have reported the discovery of stars (e .g.. Majewski et all 
2003; Martmez-Delgado et al. 2004; B elokurov et'aiTl2006l; 
I Yannv et al.. .2009b; Correnti et al.. ,20101 ; a compre hensive 
summ ary of the earlier detections appears in Majewsk i et al.l 
120031) or star clusters (e.g.. Pal 12: iDinescu et al. 2000; Whit- 
ing 1 : Carraro et al. 200"^ many clusters: Bellazzini et aH 
2003; a summary of Sgr clusters appears in "Law & Majewski 
i2010b.) plausibly associated with debris from Sgr, either 
trailing or leading it along its orbit. Line-of-sight ve- 

* Though we note that th ere is now evidence for ext ended tidal debris pop- 
ulatio ns around the Carina (Muiioz et al.'2006bV2008) and Leo I (So hn et al] 
|2007|) dSphs. Also, some debate still exists over whether the HI Magellanic 
Stream derives from tidal stripping of Small or Large Magellanic Cloud gas 
versus from ram pressure stripping. 
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Figure 1. Left panel: Distiibution of Kapteyn's Selected Areas (in equatorial coordinates) for which we have derived proper inotions, shown in an Aitoff 
projection. Solid points are those fields for which we have additional deep, 4-meter plates (see text). The current orbital plane of Sagittarius is overlaid as a solid 
blue line, and the shaded (light green) areas represent the sky coverage of SDSS (as of DR5). Regions containing stellar overdensities suggested in the literature 
to be part of the "Monoceros ring" are denoted by the blue hatched areas on either side of the disk. Right panel: Spatial distribution of the Kapteyn Selected 
Areas used in this study overlaid on the predicted distribution of Sagittarius tidal debris from the best-fit triaxial halo model of Law & Majewski (2010a). Gold 
colored points represent debris stripped from the Sgr progenitor on the past two perigalactic passages (0-1.3 Gyr ago), and magenta points the previous two 
passages (1.3-3.2 Gyr ago). Note that all of the fields are sainpling predominantly debris stripped on the saine orbital passage (i.e., the gold points), with only SA 
71 slightly sampling earlier-stripped (magenta) debris. 



locities (i.e., radial velocities, or RVs) of Sgr mem- 
bers have been determined at a few positions along the 
stream (e.g., Dohm-Palmer et al. 2001; Majewski et al. 200^ 
[Monaco et al.ll2007l) . and, along with the spatial distribution 
of these stars, provide constraints on models of the Sgr-Milky 
Way interaction (e.g., [Johnston et al. 1995; Helmi & Whi3 
2001'; 'Ibat a et al.ll2001l: lHelmiti2004i ,Martmez-De lgado et af 
2004). A comprehensive effort at modeling the Sgr disruption 
constrained by all observations available after about a decade 
of study was done bv lLaw eFal. (2005). who were able to re- 
produce most extant data, but were unable to completely rec- 
oncile the apparent need for a prolate MW halo potential to 
produce the leading arm radial velocities on the one hand, 
and an oblate halo to match the positions of leading debris 
on the other This cont radiction has apparently been recently 
resolved by iLaw et aLl (|2009), who propose that the Milky 
Way might have a triaxial halo; a comp rehensive A^-body 
model based on the best-fitting triaxial halo (iLaw & Maiewskil 
l2010al hereafter LMIO) reasonably matches nearly all exist- 
ing constraints (spatial and kinematical) of Sgr tidal debrisQ 
The Sagittarius dwarf and its tidal debris are thus proving to 
be an excellent laboratory for studying both the dynamics of 
tidally disrupting dwarf galaxies and star stream formation, 
as well as the shape and strength of the Galactic gravitational 
p otential that is the cause of this disruption. It is this model 
oflLMlOl to which we shall compare our data throughout this 
work. 

1.1. Our Proper Motion Survey 

' Further co mplications have arisen du e to an a pparent bifurcation of the 
leading stream fflelokurov et al.'200ff); the'LMlO' model was not designed to 
address this issue. Several attempts to explain this detail invoke overlapping 
debris from multiple orbital wraps I Fellhau er et al. 2006; though Yan nv et alj 
1200 9b find siinilar stellar populations in both arms, likely ruling out this 
scenario) or a disk-galaxy progenitor for Sgr (Peiiarrubia et al. 2010j; but cf. 
[Epkas et al. 2010). The highly-elliptical shape of the Sgr dwarf has recently 
been reproduced by (Lokas et al. 2010), who inodel Sgr as a disk galaxy em- 
bedded in an extended dark halo. Tidal stirring transforms the initially disk 
dwaii' into an extended elliptical shape over two pericentric passages; the ro- 
tation of the progenitor may also explain the bifurcation of the Sgr leading 
arm. 



To date, no systematic survey has addressed the tangen- 
tial velocities (derived from proper motions) of the iden- 
tified major Galactic tidal streams. Only a few studies 
(e.g., Qp inescu et al. 2002; Casetti-Dines cu et al.l l2008. 2009t 
iCarlin et al.. .2010' ,Koposov et al. .2010.) have published any 
proper motion results for major Galactic substructures, and 
typically not at a level of precision that is useful for constrain- 
ing dynamical models of tidal stream production and evolu- 
tion. In an effort to detect and characterize halo substruc- 
tures, we have been working on a project to obtain full phase- 
space information (positions and full 3-D space motions) 
for individual stars in Kap teyn's Selected Areas (SAs; see 
ICasetti-Dinescu et al.ll200d for an overview of this project). 
The sky positions of the Selected Areas were chosen by Ja- 
cobus Kapteyn ( 1906) to provide evenly spaced coverage for 
a systematic exploration of Milky Way structure. We have at- 
tempted to carry on at least part of this legacy by taking advan- 
tage of Mt. Wilson 60-inch telescope phot ographic plate ma- 
terial taken by Kapteyn and collaborators dSeares et al.lll930l) 
for their survey to make up the first-epoch data of our survey 
(in particular, near-equatorial fields at (5 = 0°, +15°, and-15°). 
The distribution on the sky of those SAs that make up our sur- 
vey is shown in an Aitoff projection in Figure [T] Some of the 
equatorial SA fields lie along the orbit of the Sagittarius dwarf 
galaxy (which is approximated by the blue, solid curve in Fig- 
ure [Til - it is a subset of these fields (in particular, six fields 
along the trailing tidal tail; see the right panel of Figure[T]i that 
are the focus of the present work. 

1.2. Constraints on the Local Standard of Rest Velocity 

The positions of tidal debris that have been found over a 
large stretch of the Sagittarius orbit place fairly strong con- 
straints on the three-dimensional motions of the Sgr dwarf. 
It is, however, important to confirm and refine the models by 
measuring space velocities of stream stars (especially proper 
motions, which are difficult to measure for stars in distant 
Galactic substructures). In the case of the Sgr trailing tail, 
however, proper motions measured for debris stars are also 
rather sensitive to the Sun's motion through the Galaxy. This 
arises because much of the Sagittarius trailing tidal tail is 
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positioned at a roughly constant distance below the Galactic 
plane, with the Sgr orbital pla ne nearly coiri c ident al with the 
Galactic Xqc - Zqc planeQ Maiewski et alj (12006 , hereafter 
"MLPP") noted that because of this orientation, longitudinal 
proper motions of Sgr trailing debris located sufficiently far 
away from the South Galactic Pole contain virtually no con- 
tribution from Sagittarius motions, and almost entirely reflect 
the solar motion. 

Efforts to measure fundamental dynamical properties of the 
Milky Way, such as its rotation curve, Q{R), are compli- 
cated by our Sun's own (poorly known) motion within the 
Galaxy. Measurements of 9lsr, the Galactic rotation speed 
at the solar circle (the Local Standard of Rest, "LSR"), vary 
by 25%, despite many efforts at its determination. The value 
adopted by the lAU in 1985 of Olsr = 220 km s"' (see 
iKerr & Lvnden -Bell 1986) has long represented a reasonable 
approximation to existing measurements (note, however, that 
prior to the 1985 lAU adoption of 8lsr = 220 km s"', the 
1964 lAU general assembly adopted 250 km s~'; see a listing 
of pr e- 1985 measurements of 9lsr in iKerr & Lynden-Belll 
Il986h . Constraints taking into account the ellipticity of the 
disk have suggested the LSR velocity could be as low as ^ 1 80 
km s~' (iKuiiken & Tremainelll994i) . A similarly low value 
of 184 ± 8 km s"' was found bv lOUing & Merrifieldl (119981) , 
who modified previous methods of determining the Oort con- 
stants by including radial variations of gas density in their 
mass modeling of the Galactic rotation cu rve. Proper motions 
of Ga lactic Cepheids from Hipparcos (iFeast & Whitelocfl 
[T997h yield a result of Olsr = (217.5 ± 7.0)(/Jo/8) km s"' 
(where Rq is the distance from the Sun to the Galactic center; 
the lAU adopted value is Rq = 8.5 kpc), in line with the lAU 
standard. Using re -reduced Hipparcos data dvan LeeuwenI 
l2007h with improved systematic errors, lYuan et all (120081) 
found Olsr = (243 ± 9)(Ro/8) km s"' based on thin-disk O- 
B5 stars. Estimates based on absolute PMs of Galactic bulge 
stars in the field of view of globular cluster M4 using the 
Hubble Space Telesc ope (HST) yield (202.4 ± 20.8)(/?()/8) 
kms"' dKalirai et al.l l2004) and (220.8 ± 13.6)(/;o/8) km s"' 
(iBedin et al.l l2003l) ~(with both studies using data from the 
same HST observations). Long-term VLBA monitoring of Sgr 
A*, the radio source at the Galactic center, led to a proper 
motion of Sgr A* from which iReid & Brunthaleii (120041) re- 
vise d the LSR velocit y upward to (235.6 ± 1.2)(7;o/8) km 
s"' . iGhez et aP (120081) combined stellar kinematics near the 
Galactic center with the proper motion of Sgr A* to derive 
(229± 18)(/?()/8.4) km s"'. More recently, Olsr has been 
suggested to be even higher, (254± 16)(/?o/8.4) km s"', based 
on trigo nometric parallaxes of Galactic star-forming regions 
(iReid e t al. 2009). Reanalysis of these same data, includ- 
ing the Sgr A* proper motion, by Bow et al. ( 2009) found 
a similar (244 ± 13) km s~'. Ko posov et a l. (2010) provided 
constraints on the MW halo potential by analysing the GD-1 
dGrillmair & Dionatosll200 6b stellar stream, combinin g SPSS 
photometry, USNO-Bh-SDSS proper motions (see Mun n et aTl 
l2004il2008h . and spectroscopy to obtain 6-D phase-space data 

^ Throughout this paper, when we refer to Galactic Cailesian (X,Y,Z)qc 
coordinates, we are specifically referring to a right-handed Cartesian frame 
centered on the Galactic center, with Xqc positive in the direction from the 
Sun to the Galactic center, Yac in the direction of the Sun's motion through 
the Galaxy, and Zqc upward out of the plane. Assuming the Sun is at = 8 
kpc from the Galactic center, this places the Sun at (X, Y, Z)gc=(-8 0,0,0) kpc. 
The con'esponding velocity components will be denoted (U ,V,W)ac^ where 
the "GC" denotes velocities relative to the Galactic rest frame. 



over a large stretch of the stream, which they used to estimate 
Slsr = (224 ± 13)(/?o/8.4) km s"' (though this result is made 
somewhat more uncertain due to a systematic dependence on 
the flattening of the disk-i-halo potential). Finally, a combined 
estimate including many of the above results as priors finds a 
value of (236 ± 1 l)(7;o/8.2) km s"' dBovv et al.ll2009l) . Most 
of the estimates discussed here rely on the Oort constants, and 
thus are dependent on our incomplete knowledge of Rq. De- 
spite numerous attempts at determining the circular velocity 
at the solar circle, this constant remains poorly constrained. It 
is clear that independent methods would be valuable to obtain 
alternative estimates of Glsr- 

Here we use a new, independent method for ascertaining 
©LSR that has the advantage over most of the previously men- 
tioned methods in that the results have virtually complete 
decoupling from an assumed value of Rq. As discussed in 
MLPP, the trailing arm of the Sagittarius tidal stellar stream 
is ideally placed to serve as an absolute velocity reference for 
the LSR. With an orbital pole of (lp,bp) = (274,-14)°, Sgr 
is almost on a polar orbit, and the line of nodes of the inter- 
section of the Galactic midplane and the Sgr debris plane is 
almost coincident with the Galactic Xqc axis (the axis con- 
taining the Sun and Galactic center). This is illustrated in 
Figure [2l which shows the projection of Sgr debris from the 
ILMIOI model onto the Galactic Xgc-Zgc, Ycc-Zcc, and 
Xgq-Ygc planes. In the upper left panel (the Xgc-Zcc plane), 
the Sgr orbital plane is nearly face-on, while in the other two 
panels, few Sgr debris points are seen more than ~ 5 kpc 
on either side of the Xcc-Zgc plane (i.e., jlGcl ^ 5 kpc for 
nearly all Sgr debris). The motions of Sgr sta rs within its (vir- 
tually non-precessing; iJohnston et al.ll2005h debris plane, as 
observed from the LSR, are therefore almost entirely in the 
Galactic U and W velocity components (i.e., in the Xqc-Zcc 
plane), whereas V motions of Sgr tidal tail stars almost en- 
tirely reflect solar motion — i.e., 9lsr (plus the Sun's pecu- 
liar motion i n V, established to be in the range ^ -1-5 to -1-12 
km s"'; e.g.. iDehnen & Binneylll99 8). The Sgr trailing tail 
is positioned fairly equidistantly from the Galactic disk for 
a substantial fraction of its stretch across the Southern MW 
hemisphere (iMajewski et al.ll2b03l) . This band of stars arcing 
almost directly "beneath" us within the Xcc-Zcc plane (see 
the upper panel of Figure|2|l provides a remarkable, stationary 
zero-point reference against which to make direct measure- 
ment of the solar motion almost completely independent of 
the Sun 's distance from the GC. 

Because of the fortuitous orientation of the Sgr debris, the 
majority of 8lsr motion (i.e., V) is seen in the proper motions 
of these stars, with the reflex solar motion almost entirely con- 
tained in the /i/ cos(b) component for Sgr trailing arm stars (at 
least for those stream stars away from the South Galactic Pole 
(SGP) coordinate "discontinuity", where the /i;cos(b) of Sgr 
stream stars switches sign). Fig. 4 of MLPP shows the essence 
of the proposed experiment via measurement of /i/ cos(b) for 
Sgr trailing arm stars, which shows a trend with debris longi- 
tude, AgU that reflects the solar motion. In the region from 
100° < A0 < 200°, /i/cos(b) is nearly constant, because the 
motion of Sgr debris contributes little to the V-component of 
velocity. Thus accurate measurement of jii cos(b) for Sgr trail- 
ing tail stars along this stretch of the stream will provide a 
means of estimating Glsr with almost no dependence on Rq. 

' Kq was defined bv lMaiewski et al]|2003l as longitude in the Sgr debris 
plane as seen from the Sun; Aq = 0° at the present Sgr position, and increases 
along the trailing tail. 
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Five of the Kapteyn fields for which we have precise (^1 
mas yr"') proper motions he squarely on the Sgr trailing arm 
in this A0 range, and o ne ot her (SA 92) is on the periphery 
of the stream. In Section l4!2l we will use the mean Sgr debris 
proper motions derived in four of these six fields to derive 
constraints on 9lsr- 

1.3. Metallicities and Detailed Abundances of the Sagittarius 

System 

IChou et aH (l2007h presented one of the first studies of 
high-resolution spectroscopic metallicities derived for Sgr de- 
bris. Their work showed that M-giants along the Sagittar- 
ius leading stream exhibit a significant metallicity gradient 
(which had previously been suggested to be present over 
smaller separatio ns from the Sgr core based on photometric 
techniques; e.g., Alard 2001; 'Mar tmez-Delgado et al.ll2004t 
[Bellazzini et al. 2006), decreasing from a mean [Fe/H] = - 
0.4 in the core to ^ -0.7 between ^ 60- 120° from the core 
(i.e., between 300 > Aq > 240°) , and to - -1.1 at > 300° 
from the main Sgr body. Such a population gradient along 
the stream likely arose due a strong metallicity gradient be- 
ing present in the dSph before its tidal disruption; thus the 
outer, more metal-poor populations were preferentially lost as 
tidal stripping progressed at earlier times relative to the more 
intermediate-age (and higher metallicity) populations remain- 
ing in the core (a mechanism for this process has been demon- 
strated in the context of an A^-body model by LMIO). In 
addition, apparently some younger populations were formed 
even after Sgr began disrupting. The existen ce of a popula - 
tion gradient has also been seen by Bellaz zini et al.l (l2006l) . 
who found that the relative numbers of blue horizontal branch 
(BHB) stars to red clump (RC, or red horizontal branch) stars 
are much higher in a leading stream field than in the Sgr core. 
Since BHB stars arise in older, more metal-poor populations 
than the RC stars, this must indicate that the stripped popula- 
tion was made up of predomina ntly older, less-enri ched stars 
than remain in the core today. iKeller et"an ( 120101) extended 
the search for chemical evolutionary signatures to the trail- 
ing tail of Sgr, observing a ha ndful of stars se l ected from the 
2MASS M-giant catalogs of Maiewski et al.1 (I2003h at high 
resolution in each of two fields at distances of 66° and 132° 
from the core. Keller et al. combined the mean metallicities in 
these two fields wit h the [Fe/H] = -0.4 result for the Sgr core 
from [Monaco et all (j2005), and derived a metallicity fit as a 
function of Aq of A[Fe/H] = (-2.4±0.3) x 10"^ dex degree"'. 
This trend (seen in their Figure 4) also pass es through the 
mean metallicity of [Fe/H] k, -0.6 derived bv iMonaco et al.l 
(l2007h in a narrow region centered at Kq = 100°. 

For consistency, all of the data i ncluded in the Keller et aD 
(|2010) study (including those from lMonaco et al.lr2005L 120071 
and Chou et al. 2010) were derived from M giants, which 
are, however, biased toward metal-rich, and therefore rela- 
tively younger, stars. In the current study, we explore the 
metallicity in fields between 75 < A0 < 130° from the Sgr 
core along the trailing tail using predominantly main sequence 
stars. Such stars near the main sequence turnoff are much 
less prone to metallicity biases than M giants, because MSTO 
stars are present in all stellar populations. An additional ad- 
vantage of focusing on MSTO stars is that the number density 
of turnoff stars is much higher than both young, M giant trac- 
ers and older horizontal-branch stars; this provides us a much 
larger sample with which to characterize the Sgr trailing tail 
metallicity. Older trailing debris populations have recently 



been studied by ISesar et all (120101) . who used SDSS Stripe 
82 data to develop a new technique for estimating metallic- 
ity from photometric data where both RR Lyrae variables and 
main-sequence stars from the same structure can be identified. 
Their work found a constant [Fe/H] = -1.20±0.1 for Sgr de- 
bris along much of the same region of the trailing tail we are 
studying. In Section 5 we explore the Sgr trailing tail metal- 
licity based on our samples of predominantly MSTO stars. 

An important diagnostic of the star-formation timescale in 
a system is the a-element abundance; the a elements (e.g., 
Mg, Ca, Ti) are produced mainly in Type II supernovae (SNe), 
which are the evolutionary endpoints of massive stars that 
dominate the chemical evolution at early times. Once Type 
la SNe begin to occur, the [a/Fe] ratio will decrease, be- 
cause a-elements are less effectively produced by these su- 
pernova progenitors, while the overall metallicity, [Fe/H], will 
continue to increase. This produces a "knee" in the [a/Fe] 
vs. [Fe/H] diagram, which acts essentially as a chronome- 
ter for a given system, since the [Fe/H] of the knee indi- 
cates the transition from SNII-dominated evolution to SNIa 
contributions. This phenomenon has been seen in a num- 
ber of dSph systems, which typically show lower [a/Fe] at 
a given [Fe/H] than Galactic populations beca use of a slower 
enrichment (e.g., Shetrone et al. 2001, 2003; Venn et al.l2004t 
iTolstovet al. 1 12003; Geisle r et al. .2005; see also a recent re- 
view bylTolstoy et al. 2009^ 

However, at the lowest metal- 
licities, the [a/Fe] of dSphs m ore closely res e mble those 
of the MW halo. Studies by iSbordone etall (I2007h and 
IMonaco et al.l (12005) found the same underabundance of a- 
elements relative to the Milky Way for the core of the Sagittar- 
ius dSph. This finding has been extended into the Sgr trailing 
stream by [M onaco et al.' (2007), and into the leading arm by 
IChou et ail (2010), with both studies using M-giants from the 
catalog of Maie wski et al.i (l2003i) . However, M giants are bi- 
ased to higher metallicity and more recently star-forming pop- 
ulation(s) of Sgr, so a natural next step in understanding the 
evolution of the original, pre-disruption Sgr dSph is to derive 
detailed abundances (especially for s-process and a-elements) 
for a significant sample of the more metal-poor, older stars 
populating the core or, more accessibly, in Sgr's more nearby 
streams. In Section 5 we present relative Mg abundances de- 
rived from our spectra. We show that the majority of con- 
firmed old, metal-poor Sgr stream members appear to have 
distinct Mg abundances from those of the Milky Way stellar 
populations along the lines of sight probed. 

1.4. Goals of This Paper 

Here, we present data in six of the Kap teyn's Selected Areas 
from our deep proper- motion survey ( Casetti-Dinescu et al.l 
120061) . In these six fields intersecting the trailing tidal tail 
of the Sgr system, we have augmented our proper-motion 
catalogs with follow-up spectroscopy. Sgr debris has been 
identified from among the stars with measured radial veloc- 
ities, and these Sgr candidates are used to derive the mean 
three-dimensional kinematics and chemistry of the Sgr trail- 
ing stream. 

In Section 12.21 we briefly introduce the proper mo- 
tion sur vey (a m ore detailed discussion of the survey ap- 
pears in iCasetti -Dinescu et al. 2006), and discuss in depth 
the spectroscopic observations with the WIYN-nHydra and 
MMT-nHectospec multifiber instruments that yielded a total 
of > 1500 radial velocities among proper motio n-selected 
stars within the six fields of view. Section 12.31 will detail 
the final selection of candidate Sgr debris in each of our 
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fields based on RV, proper motion, and color-magnitude se- 
lection. Section[3]presents maximum-likelihood estimates of 
the most precise proper motions (^1-2 mas yr"' per star, 
or ~ 0.2-0.7 mas yr"' mean for each field) yet measured 
for Sagittarius debris. These measured k inematics are com- 
pared to the models of iLaw & MajewskH (l2010ai) . and found 
to agree rather well with the predictions for Sgr debris mo- 
tions. However, we follow in Section |4] with an analysis of 
the residual disagreement between our measurements and the 
models, or more accurately, we use the discrepancy to re- 
assess the magnitude of the solar reflex motion, which is the 
dominant contributor to the proper motions in the direction 
of Galactic longitude. We show that our proper motion data 
(specifically, fiicosb) are inconsistent with the standard lAU 
value of 220 km s"' for the Local Standard of Rest motion at 
the ^1-2(7 level and favor a significantly higher value, con- 
sistent with several of the most recent 0lsr studies using ra- 
dio techniques. In Section|5]we apply a software pipeline de- 
signed to derive stellar abundances from low-resolution spec- 
tra to the numerous spectra we have obtained for this project. 
While the metallicities show a hint of a gradient among the 
metal-poor stars in our study consistent with previous work, 
we cannot rule out a constant [Fe/H] over the range of stream 
longitude covered. We also examine the relative magnesium 
and iron index strengths for information on a-abundance pat- 
terns of Sgr debris. We find Mg abundances of Sgr members 
are typically lower at a given [Fe/H] than field stars, consistent 
with the behavior seen in most MW dSphs. Finally, Section|6] 
concludes with a brief summary of our work, and future av- 
enues these data can be used to explore. 

2. THE DATA 

2.1. Field Locations 

The data discussed here are part of our ongoing deep 
proper-motion surve y in a subset of Kapteyn's S elected Ar- 
eas (iMaiewskill 19921: ICasetti-Dinescu et al.ll2006l) at declina- 
tions o f ±15° and 0°. The survey as designed by iKaptevnl 
( 119061) consists of ~ 1° fields evenly spaced at '-^ 15° in- 
tervals along strips of c onstant declination (see Fig. 1 in 
ICasetti-Dinescu et al .''2006). A handful of our near-equatorial 
survey fields (see the left panel of Figure [U fall on or 
near the location of Sgr trailing tidal debris as mapped by 
iMaiewski et~aLl (12003 ) using M-giants from 2MASS. The lo- 
cation of our fields relative to the models (constrained by 
the Majewski et al. data, among others) of Law & M ajewskil 
(pOlOa) can be seen in the right panel of Figure[T] which sug- 
gests that we can expect a significant contribution from Sgr 
de bris to the stellar pop u lations alon g these hnes of sight. 
In ICasetti-Dinescu et al.l (120061) and ICasetti-Dinescu et al] 
(120081) we showed that faint (V or ^ > 20), blue (B-V < 0.8 
OT g-r < 0.6) overdensities (at colors and magnitudes con- 
sistent with the expected Sgr main sequence turnoff) in the 
color-magnitude diagrams of those SA fields intersecting the 
Sgr orbital path also show clumping in the distributions of 
their proper motions; these excesses and their clumping in 
proper motion space suggest that a distant, common-motion 
population, likely to be Sgr tidal debris, is present in these 
fields. In this work, we focus on six fields: SAs 1 16, 1 17, 92, 
93, 94, and 71 (listed in order of increasing A0). Coordinates 
for these fields are given in Tablefl] which includes equatorial 
and Galactic positions as well as the longitude and latitude in 
the Sagittarius coordinate system. In this study, we are dis- 
cussing fields that are ~ 74° - 128° from the core of the Sgr 
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Figure 2. Sagittarius d ebris from the best-fit triaxial halo model of 
'Law & MaiewskJ '2010a'), shown in Galactic (X,Y,Z)qc (right-handed) co- 
ordinates. Colors represent debris stripped on successive orbits, as in Fig- 
ure[T] with 2 passages of additional (earlier) debris included as cyan points. 
The three panels represent the projection of Sgr debris onto the Galactic 
^GC ~Zgc, Yqc -Zqc, and Xqc -Yqc planes. In the upper left panel (the 
^GC ~Zgc plane) the Sgr orbital plane is nearly face-on; open black squares 
in this panel denote the positions of the Kapteyn Selected Ai'eas in this 
study along the Sgr trailing tail. The Sun is represented by the circle at 
(Xqc, Zqc )=(-8. 0,0.0) kpc, with the Sgr core (black dots) beyond the Galactic 
center as viewed from our position, and sHghtly below the plane. The upper 
right and lower left panels (i.e., the Iqc— Zqc and Xqc— ^gc planes) illustrate 
the near-coincidence of the Sgr orbital plane with the Galactic Xqq - Zqc 
plane. Note that very few Sgr debris points make excursions of more than 
~ 5— 10 kpc to either side of the Galactic Xqq —Zqc plane (in the Yqq direc- 
tion). 

dSph, along its trailing stream. 

The positions of the SAs in this study in the Galactic Carte- 
sian Xcc - -Zgc plane are shown in the upper left panel of 
Figure |2] overlaid atop simulated Sgr debris from the ILMToI 
model. Of course, to place points on this figure for the SAs 
requires an estimate of heliocentric distance. Where needed 
throughout this work, we use mean distances to Sgr debris in 
each SA field estimated from the LMIO model debris along 
corresponding lines of sight. We have chosen to do this rather 
than measure Sgr debris distances because (a) our data in most 
fields don't reach much fainter than the main sequence turnoff 
of Sgr, and (b) the expected line-of-sight depth of the Sgr 
stream in this portion of the trailing tail is 10 kpc, which 
"smears" the main sequence out by as much as ^ 0.5 magni- 
tudes. Both of these factors make isochrone fitting to derive 
distances rather unconstrained, so for all analysis requiring a 
dista nce esti mate we adopt the mean model debris distances 
from ILMIOI at each position, with the line-of-sight depth of 
the stream in each field defining the uncertainty in the dis- 
tance. These values are given in Table [3]below. 

2.2. Photometry and Proper Motions 

For those fields (SAs 92, 93, 94, and most of SA 116) 
that lie within the Sloan Digital Sky Survey (SDSS) foot- 
print, we have used photometry (shown in Figure O from 
SDSS Data Release 7 (DR7; Abazajian et al. 2009) for our 
analyses. The remaining photometric data for this survey 
(for SAs 71 and 117) are photographic and derived from 
the late epoch (du Pont 2.5-m telescope) plates from which 
the proper motions were measured. Calibration of the pho- 
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Table 1 

Kapteyn's Selected Areas in This Study 



SA 


RA 


Dec 


I 


h 


Aq" 


Bq 


E(B-Vf 




(J2000.0) 


(J2000.0) 


(degrees) 


(degrees) 


(degrees) 


(degrees) 




71 


03:17:11.5 


15:24:57.6 


167.1 


-34.7 


128.2 


-5.6 


0.19 


94 


02:55:58.1 


00:30:03.6 


175.3 


-49.3 


116.3 


4.8 


0.09 


93 


01:54:52.1 


00:46:40.8 


154.2 


-58.2 


103.2 


-3.2 


0.03 


92 


00:55:03.8 


00:47:13.2 


124.9 


-62.1 


90.1 


-10.6 


0.03 


117 


01:17:04.1 


-14:11:13.2 


149.0 


-75.7 


87.6 


5.1 


0.02 


116 


00:18:08.4 


-14:19:19.2 


90.1 


-75.0 


74.9 


-1.4 


0.02 



"Coordinates in the Sagittarius system as detined bv lMaiewski et al.l <2003l) . Aq and Bq are analogous to Galactic longitude and latitude, but rotated such that 
the Sgr core defines the center of the system (i.e., Aq,B0 = 0°,0°), with Aq increasing along the trailing tidal tail. The fields in this study sample the trailing 
tail between 74—128° from the Sgr core. 

''Interstellar reddening value estimated from the maps of lSchlegel et aljj 19981) . 
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Figure 3. Left: SDSS color-magnitude diagrams (CMDs) of all stars with measured proper motions in the four fields of our survey that overlap the SDSS 
footprint. Each of these CMDs shows the blue (g - r ~ 0.5) swath of stars at the bright end made up of primarily Milky Way thin- and thick-disk MSTO stars. 
Below this feature at similar blue, (g — r< 0.6) colors, but at fainter {g > 19.5) magnitudes in each field is an apparent overdensity likely made up of Sagittarius 
main-sequence stars. Note that the much deeper proper motion catalog of SA 94 samples much more of the Sgr main sequence than those catalogs for the other 
fields. Right: Proper motion vector point diagrams (VPD) , separated into a blue (0.2 < g—r < 0.8) sample (middle column) and a red (1.2 < g-r < 1.7) subset 
(right). Red stars from the prominent feature visible at red colors (g-r > 1.2) in each CMD should be primarily nearby Galactic M-dwarfs. The blue stars 
contain main sequence turnoff stars of MW populations, as well as candidate Sgr MSTO stars. The proper motions of many stars are tightly clumped in the blue 
samples, suggesting that a distant, common-motion stellar population (i.e., Sgr debris) may be present among these stars. 



tographic magnitudes in the blue (IIIa-J-i-GG385) and vi- 
sual (IIIa-F-i-GG495) passbands onto the standard Johnson- 
Cousins system was achieved using CCD photometry taken 
in 1997-1998 with the Swope 1-m at Las Campanas Obser- 
vatory. However, that UBV CCD photometry only covers a 
small portion 20 - 30%) of each field, and is shallower 
than the magnitude limit of the photographic plates by 1- 
1.5 magnitudes for red stars, and 0.3 magnitudes for blue 



stars, yielding poor photometric calibration (for details, see 
Casetti-Dinescu et al. 2006). Thus, for SAs 71 and 1 17, those 
aspects of our analysis that rely on the photometry (e.g., spec- 
troscopic target selection, photometric parallax distance esti- 
mates, and calibration of color-dependent systematics in the 
proper motions) are subject to the uncertainties in calibra- 
tion of the photographic photometry propagated by up to 0.05 
magnitudes in the B and V photometry. 



Kinematics of Sgr Tidal Debris in Kapteyn's Selected Areas 



7 



Details of the p r oper motion reductions appear in 
ICasetti-Dinescu et alj (l2006h . so here we provide only an 
overview. For all of the near-equatorial (-15° < S < 15°) Se- 
lected Areas in our study, proper motions are derived from 
plates taken with the Mt. Wilson 60-inch between 1909-1912, 
combined with deliberately matched plates (in approximate 
area and plate scale) taken by S. Majewski with the 2.5-m Las 
Campanas du Pont telescope between 1996-1998. All plates 
were digitized with the Yale PDS microdensitometer Most 
background QSOs and galaxies are near the limiting magni- 
tude of the proper motion catalogs derived from solely the 
Mt. Wilson and du Pont plates; therefore if we used only these 
data, the correction to an absolute proper motion frame would 
be determined by only a handful of poorly-measured faint ob- 
jects. To extend the proper motion limiting magnitude beyond 
the limit imposed by the Mt. Wilson plates, we augmented 
the Kapteyn survey data with measurements of plates from 
the first Palomar Observatory Sky Survey (POSS-I), which 
were taken in the early 1950s. These plates, while of much 
coarser plate scale (67 ."2 mm"'; compare to 10"92 mm"' for 
the du Pont, and 27" 12 mm"' for the 60"), are deeper than 
the 60" plates, and provide a ^ 40-year baseline with both 
the Mt. Wilson plates and the du Pont plates, allowing us to 
extend the limiting magnitude of the survey (at least for the 
POSS-I/du Pont proper motion baseline) to V > 20.5. 

In addition, for two of the S A fields in this study (S As 7 1 
and 94) we have included plates from the Kitt Peak National 
Observatory Mayall 4-meter telescope, taken at prime focus in 
the mid-1970s by A. Sandage and in the mid-1990s by S. Ma- 
jewski. This additional high-quality plate material extends the 
limiting magnitude in SAs 71 and 94 to ~ 50% completeness 
at SDSS r magnitudes of ^ 22 (compared to ^ 50% complete- 
ness at r ~ 19.5 in fields without 4-meter plates; the difference 
in depth between SA 94 and the other fields can be seen in 
Figure O. The precision of the proper motion measurements 
for the "deeper" fields (i.e., those with 4-meter plates) is im- 
proved by roughly a factor of 2 at the same magnitude over 
those with no 4-m plate material; this arises because of both 
more plate material with fine (18. "6 mm"') plate scale at inter- 
mediate epochs and the increased depth provided by the 4-m 
plates. 

2.3. Radial Velocities 

The survey fields in which we focus this Sg r study fall on or 
near the portion of the trailing stream in which lMaiewski et al.l 
(|2b04) and Monaco et al. (2007) have identified a clear Sgr ra- 
dial velocity signature. These fields can be seen relative to the 
orbital path of the Sgr dSph in the left panel of Figure [T] and 
with respect to the expected location of Sg r trailing tidal de- 
bris acc ording to the best-fitting models of iLaw & Majewski' 
( |2010a|) i n the right pan el of Figure[T] We have already shown 
evidence ( Casetti-Dinescu et al.l200 6, 2008) that the overden- 
sities of faint, blue stars that are tightly clumped in proper 
motions in a few of these fields are likely made up of Sgr 
debris. It was these apparent overdensities that guided our 
target selection for spectroscopic follow-up. We began with 
spectroscopy from the Hydra multifiber spectrograph on the 
WIYN 3.5-m telescope; in most of the shallower fields of this 
survey, moderate-resolution spectra can be obtained with this 
instrument in a reasonable amount of observing time. For the 
deep fields (and some of the shallower fields as well), we used 
another multi-object spectrograph, the Hectospec instrument 
on the MMT 6.5-m, which allowed us to observe > 200 Sgr 



stream candidates simultaneously per setup down to faint (g 
or y > 21.5) magnitudes. We describe the observations and 
data reduction for each instrument separately below. 

2.3.1. Sample Selection 

Targets for spectroscopic follow-up were selected to lie 
within the locus of the suspected Sgr main sequence turnoff 
(MSTO) at faint (g or V > 19. 5) magnitudes and blue (g-r 
or B-V < 0.8) colors. ICasetti -Dinescu et al. (2006) showed 
that the proper motions of these Sgr MSTO candidates clump 
more tightly than those of the predominantly nearby M- 
dwarfs at red colors. The tight clumping in the proper motion 
vector point diagram ( VPD) of stars in the MSTO feature was 
used to define a selection box in proper motion space which 
should contain any Sgr debris that is present in each field, and 
eliminate a good fraction of unrelated stars of similar color 
and magnitude. Care was taken not to be too stringent with 
either the proper motion or photometric criteria, to preserve 
as many potential Sgr stars in the wings of the distributions 
as possible. Because the quality and depth of the photometry 
and proper motions varies between fields, different candidate 
selection criteria were adopted for each field. 

For WIYN-nHydra observations, only stars brighter than 
20th magnitude (either V or g, depending on whether a given 
field had SDSS photometry) were included in the multifiber 
setups, because fainter stars than this require rather long ex- 
posures with a 3.5-meter telescope to achieve adequate signal- 
to-noise for radial velocity measurement. After all available 
fibers were filled with MSTO candidates the remaining fibers 
were assigned to targets at relatively bright magnitudes (< 18) 
that also reside within the VPD selection criteria (i.e., po- 
tential Sgr RGB stars). Stars fainter than 20th magnitude 
were targeted with the 6.5-meter MMT telescope, which eas- 
ily reaches Sgr MSTO candidates in these fields. 

2.3.2. WIYN+Hydm Observations 

Spectroscopic data were obtained during a total of eight ob- 
serving runs with the WIYN 3.5-m telescop43 between De- 
cember 2002 and November 2008. We used the Hydra multi- 
fiber spectrograph in two different setups. The first one (Dec. 
2002, Nov. 2003 observing runs) used the 800@30.9 grating 
with the red fiber cables and an order centered in the neigh- 
borhood of the Mg triplet (5170 A) and covering about 980 
A of the spectrum. This setup delivered a dispersion of 0.478 
A pix"' and a resolving power R ^ 5400 (resolution ^ 1 A). 
The second spectrograph configuration (Dec. 2005, Oct./Dec. 
2006, Oct./Dec. 2007, and Nov. 2008) used the 600 @ 10.1 
grating with the red fiber cable to yield a wavelength cover- 
age A = 4400-7200 A at a dispersion of 1.397 A pix"', for a 
spectral resolution of 3.35 A (/? ~ 1500 at A = 5200 A). This 
spectral region was selected to include the H/3, Mg triplet, 
Na D, and Ha spectral features. Typically 60-70 targets were 
placed on Hydra fibers, with the remaining 15-20 fibers placed 
on blank sky regions to allow for accurate sky subtraction. 
Each of the 2005-6 datasets was obtained in less than optimal 
conditions, including substantial scattered moonlight in Dec. 
2005 and cloudy conditions in both 2006 runs. The major- 
ity of the 2007 and 2008 data were obtained under favorable 
conditions. We further note that the Nov. 2008 observing 

The WIYN Observatory is a joint facility of the University of 
Wisconsin-Madison, Indiana University, Yale University, and the National 
Optical Astronomy Observatory. 
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Table 2 

Summary of Spectroscopic Observations 



SA Date Telescope/Instrument Exposures A'stais Mag. limit 

(seconds) 



71 


Dec 2002 


WIYN+Hydra" 


4 X 1 800 


37 






i>UV ZUUj 


vv 1 1 iNTnyuid 


1 Y 1 son A Y 1 son 

iU X loUU, H X ioUU 


74 


18 to 


_ 


Dec 2005 


WIYN+Hydra'' 


6x 1800 


27 


19.5 


- 


Oct 2006 


WIYN+Hydra*" 


6x 1800 


47 


19.0 




Uec iUUo 


WIYJN+Hydra 


4x1 oOO 


58 


19.3 




Dec 200/ 


WIYJN+Hyara 


X 2400 


53 


19.7 




Uec iOOo 


MM 1 +Hectospec 


1 o An ^ 1 or\r\ 

X ISOO, D X IsOO 


176, 193 


22.0 






... TOTAL 




503 




92 


Oct 2006 


WIYN+Hydra'' 


6x 1800 


39 


19.4'' 






SDSS 




zll 


Zl.j 






TOTAL 




254 




93 


Oct 2006 


WIYN+Hydra'' 


6x 1800 


55 


19.8'' 


- 


Oct 2007 


WIYN+Hydra'' 


6x 1800 


51 


19.9 






MlVTT-i-Hprtncnpr 
iviivi 1 T^mcuiuftucu 


d Y 1 SOD 


195 


21.5 






SDSS 




101 


21.8 






TOTAL 




292 




94 


Oct 2007 


WIYN+Hydra'' 


8x 1800 


15 


19.0'' 




Dec 2007 


WIYN+Hydra'' 


6 X 2400 


53 


19.9 




Dec 2008 


MMT+Hectospec 


8x 1800, 6 X 1800 


189, 165 


22.5 






SDSS 




88 


22.5 






... TOTAL 




432 




116 


Oct 2007 


WIYN+Hydra'' 


6x 1800 


25 


16.7'' 




Dec 2007 


WIYN+Hydra'' 


4 X 2700 


59 


17.7 




Nov 2008 


WIYN+Hydra'' 


7x 1800 


60 


19.9 






... TOTAL 




122 




117 


Oct 2007 


WIYN+Hydra'' 


8x 1800 


47 


19.6'^ 




Dec 2007 


WIYN+Hydra'' 


4 X 2400 


49 


19.7 




Dec 2008 


MMT+Hectospec 


(4x 1800)+(4x 1500) 


182 


20.5 






... TOTAL 




206 





^These WIYN+Hydra observations used the 800@30.9 grating with the red fiber cables, centered at ~ 5200 A, yielding 1.0 A resolution spectra. 

''These WIYN+Hydra observations used the 600@10. 1 grating with the red fiber cables, yielding spectra covering wavelengths from ~ 4400 — 7200 A at 

~ 3.35 A resolution. 

'^"Roughly calibrated" V magnitudes (see lCasetti-Dinescu et aU2006D . 
''SDSS g magnitude. 



run occurred after the WIYN Bench Spectrograph Upgrade, 
which included the implementation of a new collimator into 
the Bench configuration, as well as a new CCD that delivers 
greatly increased throughput. 

Table |2] summarizes the observations. Each Hydra con- 
figuration was exposed multiple times (usually in sets of 30 
min. exposures) to enable cosmic ray removal. Standard pre- 
processing of the initial two-dimensional spectra used the CC- 
DRED package in IRAFOD Frames were summed, then 1 -D 
spectra were extracted using the DOHYDRA multifiber data 
reduction utilities (also in IRAF). Dispersion solutions were 
fitted using 30-35 emission lines from CuAr arc lamp expo- 
sures taken at each Hydra configuration. On each observing 
run we targeted a few bright radial velocity standards cover- 
ing spectral types from F through early K (both dwarfs and 
giants), each through multiple fibers, to yield multiple in- 
dividual cross-correlation template spectra. These RV stan- 

" IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Reseai'ch in Astron- 
omy (AURA) under cooperative agreement with the National Science Foun- 
dation. 



dard spectra were cross-correlated against each other using 
the IRAF tool FXCOR to determine the accuracy of the ve- 
locities and remove any outliers (i.e., those stellar spectra that 
yield unreasonable cross-correlation results due to template 
mismatch or some defect, such as a poorly-removed cosmic 
ray, bad CCD column, or other unknown culprit). Measured 
velocities of the RV standards typically agreed with published 
lAU standard values to within 1-2 km s"'. Radial velocities 
for program stars were derived by cross-correlating all object 
spectra against all of the standards taken on the same observ- 
ing run. To maximize the S/N in faint, metal -poor stars, only 
^ 200 A-wide regions centered on the H/3, Mg triplet, and Ha 
absorption lines were used for cross-correlation. 

Radial velocity uncertainties were derive d using the 
IVogt et akl (1995) method, as described in i Munoz et alJ 
(2006a) and Frinchabov et al. (2006). The Tonry-Davis ra- 
tio (TDR; Tonry&DavisiJ979 ) scales with S/N, such that 
individual RV errors can be calculated directly from the TDR, 
provided you have multiple observations at varying S/N of 
some particular standard star to map the dependence. We 
have used this technique for aU datasets except those from 
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the Dec. 2006 observing run, when only a total of four RV 
standard spectra were taken. For the S A 7 1 configuration ob- 
served on this run, the RV uncertainty for each stars is derived 
as the standard deviation of the RV results from that star's 
spectrum using cross-correlation against each of the four stan- 
dards. Typical RV uncertainties for individual measurements 
for all fields were ay « 5 - 10 km s~' , with most spectra hav- 
ing S/N 15-20 per Angstrom. The uncertainty depends on 
S/N (essentially magnitude) for all stars from a single Hydra 
pointing; however, the varying exposure times between Hydra 
setups and changing observing conditions mean that ay is not 
strictly a function of magnitude in our final catalogs. 

2.3.3. MMT+Hectospec Observations 

To obtain spectra of fainter Sgr MSTO candidates, we were 
granted three nights of qu eue-scheduled NOAO observing 
time on the MMT 6.5-m[3 A total of six observing config- 
urations were obse rved with the 3 00-fiber Hectospec multi- 
fiber spectrograph (iFabricant et al.l 2005) mounted at the f/5 
focus of the MMT. Targets were selected from among the ap- 
parent stellar overdensities of blue MSTO canidates at mag- 
nitudes too faint {g > 20) to be reasonably observed with 
WIYN-nHydra, but using the same proper motion criteria used 
to choose Hydra targets. In each of these configurations, a few 
targets previously observed with Hydra were included for a 
radial velocity consistency check, and any fibers unable to be 
filled with faint stars were assigned to brighter [g < 20) Sgr 
RGB candidates. 

We used the 270 gpm grating, centered at ^6400 A, to give 
a working wavelength range of 3700-9150 A at a dispersion 
of ^^1.21 A pix"' (4.85 A resolution). This low resolution al- 
lows us to obtain adequate signal-to-noise {S/N > 10 per A) 
spectra of stars as faint as g = 22.5 to measure radial veloc- 
ities in 3-4 hours of total exposure time. We observed two 
multifiber setups each in S As 7 1 and 94, the two fields in this 
study that have 4-meter plates (and thus deeper proper-motion 
catalogs). The rest of the time was devoted to one configura- 
tion each in SAs 93 and 1 17. Each configuration had 200-230 
targets assigned (essentially as many as the fiber assignment 
software would allow), with the remainder designated as sky 
fibers. The sky fibers were chosen in areas with the nearest 
^ < 22 star at least 15" away, distributed throughout the field 
so that a number of them would fall within each of the two 
CCD chips of the Hectospec system. The number of expo- 
sures in each field, each exposure time, the number of stars 
with measured radial velocities, and the limiting magnitude 
of each spectroscopic field are given in Table |2] 

The Hectospec data were reduced using an external version 
of th e SAO "SPECROAD" reduction pipeline (Mi nketalj 
I2007h written by Juan Cabanela and called ESPECROADQ 
The pipeline automates many reduction steps, including bias- 
correction, flat-fielding, cosmic -ray rejection, fiber-to-fiber 
throughput adjustments, and sky subtraction. Wavelength cal- 
ibration was performed manually using sets of three combined 
HeNeAr calibration lamp exposures from each night of ob- 
serving. 

We derived RVs using the IRAF task FXCOR to cross- 

Observations reported here were obtained at tlie MMT Observatory, 
a joint facility of the Smithsonian Institution and the University of Ari- 
zona. MMT telescope time was granted by NOAO (proposal ID 2()08B- 
0448), through the Telescope System Instrumentation Program (TSIP). TSIP 
is funde d by the NSF. 

http : //iparrizar . mnstate ■ eciu/~ juan/research/ESPE 



correlate object spectra against fourteen template RV standard 
spectra of nine different stars ranging in spectral type from 
F through K. We first correlated the standard spectra against 
each other, and found that our measurements agree with pub- 
lished RV£3 to within 4.9 km s~' for all of these stars, with 
zero mean offset, and ctav = 3.0 km s"'. To minimize the ef- 
fects of noise for spectra with lower S/N, the cross-correlation 
was restricted to the regions around the Ha, Mg triplet, and 
H/3 fines. 

Each of the object spectra was cross-correlated against all 
14 standards, and the mean RV from each of these 14 mea- 
surements was adopted as the final result. For the queue- 
scheduled Hectospec observations, we relied on the queue to 
provide radial velocity standards. We were thus unable to ob- 
tain repeated exposures of the same RV standard star s to all ow 
us to use the "Vogt method" (as described in Section l2.3.2b to 
determine RV errors. Instead, uncertainties were estimated 
as the standard deviation of the 14 independent RV measure- 
ments thus derived, and vary (essentially as a function of S/N) 
from (Tvheiio ^ 3 km s"' at g =18.0 to ~15 km s"' at g = 21.5- 
22. 

From repeat measures of a handful of stars, including mul- 
tiple Hydra or Hectospec observations as well as many ob- 
served with both systems, we found mean systematic offsets 
of < 5 km s"' between observing runs (including both Hydra 
and Hectospec data). These offsets were applied to all RVs 
from a given run to place all measurements on the same sys- 
tem as the Dec. 2007 WIYN-nHydra velocities. 

We also note that because we selected most Hectospec tar- 
gets to be faint, blue objects with miniscule proper motions, 
a large number of obvious QSOs and AGN spectra appeared 
in our data. These were added to the samples of QSOs and 
galaxies that provided the fixed absolute proper motion frame 
in each field, improving the zero points in these fields. 

2.3.4. SDSS spectra 

We supplemented our database of radial velocities by 
matching our proper motion catalogs to the SDSS spectro- 
scopic database. The number of additional RVs contributed 
by SDSS in each field is noted in Table |2] The majority 
of SDSS stars are red, nearby M-dwarfs, so very few addi- 
tional Sagittarius candidates were contributed by the addition 
of these spectroscopic data. However, the handful of Sgr stars 
that are present, as well as any other stars in common with 
our observations, were used for a consistency check. From 
the stars in common between SDSS and our observations in 
SAs 94 and 93, we find mean offsets of < 5 km s~'. The ac- 
curacies of SDSS radial velocities ar e ~ 4 km s~' at g < 18, 
decreasing to ~ 15 km s~' at g 20 dYannv et al.ll2009ah . so 
we choose not to offset the RVs . 

3. SAGITTARIUS TIDAL DEBRIS KINEMATICS 

3.1. Radial Velocities 

Figure |4] shows all measured velocities in each of the 
six survey fields (the total number of stellar radial veloci- 
ties in each field is given in Table |2]i in the Galactocentric 
(VgseE3) frame. These consist of all RVs from WIYN-nHydra, 

From the Geneva Radial- Velocity Standard Stars at 
|htt:p : //obswww . unige ■ ch /~udry/ std/ 

Vgsr = V'helio + 9.0cosfecos/+232.0cosfosin/ + 7.0sinfo, where Vhciio is 
the measured heliocentric radial velocity. This calculation assumes a circu- 
lar velocity of 220 km s"' at the solar circle, and solar peculiar' motion of 
iW(j,A'**-)pi^-0, 12.0,7.0) km s"'. 
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Figure 4. Measured radial velocities (relative to the Galactic Standard of 
Rest) in each of the six fields, displayed in order (from top to bottom) of 
Sgr longitude (i.e., angular distance from the Sgr core), Aq. The prominent 
peak in each field at Vqsr ~ 0-100 km s"' is made up primarily of Galactic 
stellar populations. In SAs 1 17, 93, 94, and 71, an additional peak at lower 
velocities is visible. As shown in Section 3, this peak can be attributed to the 
presence of Sgr tidal debris in these fields. 

MMT+Hectospec, and where available, SDSS. Where mul- 
tiple measurements exist, the final catalog reflects the error- 
weighted mean radial velocity. In each of these fields, a broad 
peak is seen at Vgsr ^0-100 km s"' which is made up pri- 
marily of Milky Way stellar populations in these high-latitude 
fields. An additional velocity peak is evident in SAs 71, 94, 
93, and 1 17, well separated (in all fields except SA 1 17) from 
the Galactic distribution in each of these fields; it is this ad- 
ditional peak we shall show to consist of mainly Sagittarius 
trailing tidal debris. There is no readily apparent peak at 
lower Vgsr values in SAs 92 and 1 16 - this arises for differ- 
ent reasons in the two cases. SA 1 16 is the field in which we 
have the fewest measured radial velocities, and even the stars 
for which we do have data are not optimally selected to find 
Sgr debris. Because of limitations on exposure times due to 
weather, the observed stars were all at relatively bright mag- 
nitudes {g < 18), where only very few Sgr red giants should 
be found. More spectra of stars at faint > 19) magnitudes 
and blue colors will be necessary to identify Sgr debris among 
the SA 116 data. In SA 92, on the other hand, nearly twice as 
many spectra are available than in SA 116, and mostly at rela- 
tively faint magnitudes. The paucity of obvious Sgr debris in 
this field is because of the location of SA 92 on the periphery 
of the stream, where Sgr stellar densities are rather low. There 
are a handful of stars at low (Vgsr < km s~') velocities, but 
hardly enough candidates to assert that a clear Sgr presence is 
indicated in SA 92. 



To assess whether these apparent velocity overdensities are 
expected among Galactic populations in each line of sight, we 
compared the radial velocity distributions to those from the 
Besangon Galaxy model (Robin et al. 2003)0 In each SA 
field, the model query was run five times to smooth out the 
finite sampling statistics in each individual model run. The 
five catalogs were concatenated, then for each Kapteyn field 
the measured velocities were compared to the expected ra- 
dial velocities of smooth Galactic populations by scaling the 
summed Besan^on model to match the total number of stars 
in the broad peak in each RV histogram. This was done sepa- 
rately for "bright" and "faint" samples in each field, since stars 
in different magnitude ranges preferentially sample differ- 
ent Galactic populations with different velocity distributions 
(e.g., "faint" blue Galactic stars in the region of the CMD 
where Sgr MSTO stars reside will be predominantly halo 
stars, and thus have a much higher velocity dispersion than 
stars of similar color, but much brighter magnitude, where 
thin/thick disk MSTO stars predominate). 

The resulting scaled model distributions are shown as grey 
filled histograms in Figure|5]for SAs 71 and 94, with the mea- 
sured heliocentric radial velocities given as solid-lined his- 
tograms. The broad peak is reproduced well by the model 
populations, suggesting that (a) there are no large global ve- 
locity offsets present in our data, and (b) the prominent peaks 
are indeed due to foreground/background Milky Way stars. 
The additional peaks at Vheiio ~ -170 km s"' (SA 71) and 
Vheiio ~ -150 km s"' (SA 94) are clearly not due to any ex- 
pected Galactic populations along these lines of sight. 

Similar histograms are shown for SAs 93 and 117 in Fig- 
ure |6] which again clearly show that the broad, prominent 
peak in each field is made up of Galactic populations, and 
the peaks at Vheiio ~ -160 km s"' (SA 93) and Vheiio ~ -100 
km s"' (SA 117) are inconsistent with expected Milky Way 
velocities. Note that the peak in SA 117 overlaps the wings 
of the Galactic distribution, making it slightly more difficult 
to isolate bona fide Sgr members in this field on the basis of 
radial velocities alone. 

Finally, we performed the same examination in SAs 92 and 
116, with the results shown in Figure |7] Nearly all of the ve- 
locities shown in SA 92 are from the SDSS database, and are 
predominantly very red M-dwarfs. For this reason, the long 
tail of the RV distribution at negative velocities, which is due 
to thick disk and halo MSTO stars, is not well reproduced by 
our data set. There are a small number of stars at Sgr-like ve- 
locities in Figure |7] but these fall within the expected locus 
of MW stars, so we cannot definitively say that Sgr members 
are present among our SA 92 sample. In SA 116, no excess 
peak of measured RVs relative to the model predictions is ap- 
parent. This is not surprising given (a) the caveats in the first 
paragraph of this s ection regarding the data in SA 116, and 
(b) the fact that the lLaw & MaiewskH (I2010al) model predicts 
Sgr debris in this field to have RVs of -100 < Vheiio < -50 km 
s"', overlapping the wings of the Galactic distribution in this 
field. A handful of Sgr members may thus be present among 
our velocities, but they are difficult to distinguish from the 
Milky Way halo stars by their RVs. 

With RVs in hand, a next culling for Sgr stream candidates 
was obtained by simply taking all stars within a generously 
defined range around the evident associated radial velocity 
peak. Such a selection will include a few Milky Way inter- 

Model query available at |http://model.obs-besancon.ft7] 
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Figure 5. Measured heliocentric radial velocities (solid-line histograms) in SAs 71 (left) and 94 (right), divided into a bright (SA 71: V < 20.0; SA94: g < 19.5) 
sample (upper panels) and a faint (SA 71: V > 20.0; SA 94: g > 19.5) group (lower panels). The filled gray histogram in each panel is made up of Besanfon 
model points along the corresponding Une of sight in the same magnitude ranges as the data histograms, scaled to contain the same total number of stars between 
-50 < Vheiio < 100 km s"' as the observed data. In the lower panels, a Gaussian representing the best-fitting radial velocity and dispersion of Sgr candidates is 
shown (dotted curve), along with the sum of this Gaussian and the Besanfon distribution (dot-dashed histogram). These two fields include 4-meter plates in the 
proper motion deiivation, and thus contain the deepest proper motions (and the most Sgr candidates) of any fields in the survey. Very little, if any, Sgr debris is 
evident in the bright samples. 
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Figure 6. Same as Figure|5] but for SAs 93 (left) and 117 (right), and with slighdy different definitions of the bright (SA 93: g < 19.0; SA 117: V < 20.0) 
sample (upper panels) and faint (SA 93: g > 19.0; SA 117: V > 20.0) group (lower panels). In both of these fields, there is a hint of a peak at Sgr-like velocities 
in the bright samples, suggesting that a few Sgr red giants may have been identified in these fields. 



lopers, so we examined the samples thus selected to remove 
non-Sgr stars. We first removed all stars with proper motions 
|/Lt| > 10 mas yr"' in either dimension; such stars, if actually 
at the distance of the Sgr trailing tail in this region of sky 
(~25-40 kpc), would have unreasonably large (> 1000 km 
s"') tangential velocities (Vtan = 4.74c//i km s"', where d is the 
distance in kpc and /i the proper motion in mas yr"'). Faint, 
blue stars with proper motions of this magnitude must there- 
fore be nearby (foreground) MW white dwarfs or metal-poor 
subdwarfs. After removing these stars, we then examine the 
positions of all selected candidates in the color-magnitude di- 
agram. We reject faint stars that are well redward of the read- 



ily apparent Sgr main sequence, and at brighter magnitudes, 
we remove only stars at positions obviously inconsistent with 
being Sgr red giants or horizontal branch stars. In SA 93, 
a clear offset was visible between mean proper motions of 
bright (g < 19.0) Sgr candidates and fainter ones, so we chose 
to keep only candidates a.tg> 19.0, on the assumption that the 
density of stream stars should be much greater at fainter mag- 
nitudes near the lower RGB and MSTO than along the upper 
RGB. This yields fewer total Sgr candidates, but the ones that 
remain have much higher probability of being Sgr members 
than do the brighter candidates. 
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Figure 7. Same as Figure IHbutforSAs 92 (left) and 116 (right), and with slightly different definitions of the bright (SA 92: g < 19.0; SA 116: g < 18.5) sample 
(upper panels) and faint (SA 93: g > 19.0; SA 1 17: g > 18.5) group (lower panels). No excess peak is evident (relative to the Besanjon predictions) in either of 
these fields, meaning we have likely sampled very few Sgr radial velocity members along these lines of sight. 




Figure 8. SDSS g vs. g-r CMDs of SA 94; the left panel shows all objects in 
our proper motion catalog that were flagged as stars by the SDSS star/galaxy 
separator. The right panel overplots all stars observed spectroscopically as 
larger symbols: candidates within the initial RV selection are black squares 
(filled squares: final Sgr candidates; open squares: in RV selection, but re- 
moved by other criteria), open diamonds are stars with RVs outside the Sgr 
RV selection, and open circles are stars that only have RVs in SDSS (note 
that none of these ended up being selected as Sgr candidates). The final Sgr 
candidates that we selected by RV, proper motion, and color-magnitude po- 
sition (filled squares) are concentrated around a likely MSTO of Sgr debris. 
The blue curve is a Girardi et al. (2004) isochrone for a 10 Gyr population at 
[Fe/H] = -1.3 and a distance of 30 kpc. 



3.2. Selecting Final Sgr Candidates 

We now discuss how we pared down the samples of Sgr 
candidates from the initial broad RV and proper motion selec- 
tions in each field to the final, more securely-identified sam- 
ples used for analysis. For brevity, we will show detailed ex- 
amples for only two of the six fields in the study. The first 
of these is SA 94, which is the "best-case" field in our study, 
because it has deep proper motions due to the availability of 
4-meter plates in its data set. For comparison, we follow the 
SA 94 discussion with details of SA 93, which has a much 
shallower proper motion catalog than SA 94, but also has 
high-quality SDSS photometry. These two fields were cho- 
sen simply to give the reader an idea of the type and quality 
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Figure 9. Proper motion vector point diagram of SA 94, with panels and 
symbols as in Figure [s] 
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Figure 10. Reduced proper motion diagram (RPMD) for SA 94, where Hg = 
g + 5 log + with 1^1 in a rcsec yr~*. P anels and symbols are as in Figure[8] 
The blue curve is a Girardi Tt alj <2004| ) isochi'one for a 10 Gyr population at 
[Fe/H] = -1.3 and a distance of 30 kpc, with the measured mean proper motion 
of Sgr debris in SA 94 (/itot=2.34 mas yr"' ) used to convert to reduced proper 
motion. 



of the data included in this study, and the process we followed 
to select Sgr candidates. 
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3.2.1. Selecting Sgr Candidates in SA 94 

Figure [8] shows the SDSS g-r vs. g color-magnitude di- 
agram (CMD) for SA 94. The left panel shows all stars 
for which we have measured proper motions, with SDSS- 
classified galaxies removed from the sample. On the right- 
hand side, open squares depict all stars observed spectro- 
scopically to illustrate the candidate selection. The major- 
ity of spectroscopic targets in this deep proper-motion field 
were selected from the Sgr MSTO feature of faint, blue stars, 
and observed with the large-aperture MMT 6.5-m telescope. 
The remaining targets are either (a) bright stars observed with 
WlYN-nHydra, or (b) targets chosen to fill unused fibers af- 
ter all possible MSTO candidates had been assigned. The 
final culled sample of Sgr candidates (based initially on RV 
selection, with further interactive proper-motion and CMD se- 
lection performed as discussed below) is shown by the large, 
filled black squares. As expected, these concentrate at the 
Sgr MSTO locus, with a handful of brighter stars having 
properties consistent with Sgr membership as RGB or red 
clump stars. We have overlaid an isochrone from lGirardi et alJ 
(|2004) for an old (10 Gyr), metal-poor ([Fe/ H1 = -1.3)popula- 
tion at the expected distance (d = 29.5 kpc: lLaw & Majewskil 
l2010ah of Sgr debris in SA 94; the age and metallicity of this 
isochrone is chosen to m a tch th e Sgr metal-poor population 
identified by ISiegel et"an (120071) . which should be the dom- 
inant contributor to debris in this portion of the trailing tail. 
The final set of Sgr candidates in this field concentrate near 
this ridgeline; most of the scatter about the isochrone is likely 
due to the ±5 kpc line-of-sight depth of the Sgr stream as well 
as the intrinsic population dispersion in Sgr. The proper mo- 
tion vector point diagram (VPD) is displayed as Figure|9] with 
panels and symbols the same as in Figure [8] Sgr candidates 
clump more tightly in the VPD than the overall population of 
spectroscopic targets, which is consistent with the notion that 
we are indeed measuring a distant, common-motion grouping 
of stars. 

To further cull the sample, we turn to the reduced proper 
motion diagram (RPMD). Reduced proper motion, defined as 
Hg = g + 5 log 11 + 5, where g is the apparent magnitude and 
/i is the total proper motion in arcsec yr"', compresses stars 
with common tangential velocity into coherent features in the 
RPMD. Thus, a common-motion population should form a 
sequence in the RPMD, eve n if th e population has a signifi- 
cant line-of-sight depth (see iMaiew ski 1999, especially Fig- 
ure 4). In Figure[TO]we show such a diagram for SA 94, with 
the same isochrone as in Figure |8] shifted to the measured 
tangent ial v elocity of SA 94 Sgr debris (to be discussed in 
Section [33] l. After an initial calculation of the mean motion, 
candidates that were obviously inconsistent with a broadly- 
defined region (>0.3 magnitudes in Hg and/or r) about the 
ridgeline were manually removed from the sample. 

3.2.2. Selecting Sgr Candidates in SA 93 

SA 94 is one of only two fields (with SA 71) of the six 
in this study that have KPNO 4-m plates, and thus have 1-1.5 
magnitude deeper proper motions. Furthermore, of those two, 
SA 94 has higher-quality photometry (from SDSS) than the 
photographic magnitudes used for S A 7 1 . Thus S A 94 is our 
best field in terms of overall data quality. To illustrate how a 
more typical field compares to the highest-quality SA 94 data, 
we show in Figures [TT1[T2] and[T3]the same type of plots as in 
Figures [8] m and[IOl but for SA 93. This field is located in a 
high stellar density region of the stream, as is SA 94, but has 
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Figure 11. SDSS e vs. g-r CMDs of SA 93, with panels and symbols the 
same as in Figure[8| Because the proper motion catalog in this field doesn't 
reach nearly as faint stars as in SA 94, far fewer Sgr candidates have been 
identified among our radial velocities. Note also that for reasons discussed 
in the text, only stars wi th g > 19.0 were kept in the final sample. The blue 
curve is a Giraidi et al] j2004l) isochrone for a 10 Gyr population at [Fe/H] 
= -1.3 and a distance of 28 kpc. To illustrate the possible presence of very 
old, evolved stars among the Sgr candidates, we also overlay a grey dashed 
curve for the same distance, but [Fe/H] = -1.7 and an older (^^15.8 Gyr) 
population. Note that many of the open squares (RV candidates not included 
in the final sample) lie along the asymptotic giant branch of this old, metal- 
poor population. 
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Figure 12. Proper motion vector point diagram of SA 93, with panels and 
symbols as in Figure [8] 
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Figure 13. Reduced proper motion diagram (RPMD) for SA 93, where Hg = 
g + 5 log ^ + 5, with fi in arcsec yr"' . Panels and symbols are as in Figure[8] 
The blue curve is a Girardi et al. (2004) isochrone for a 10 Gyr population at 
[Fe/H] = -1.3 and a distance of 28 kpc, with the measured mean proper motion 
of Sgr debris in SA 93 (/Jtot=2.69 mas yr"' ) used to convert to reduced proper 
motion. 

shallower proper motion data (see Fig.O, providing far fewer 
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Figure 14. Proper motions of only the final Sgr candidates in SAs 94 (left 
panel) and 93 (right panel), with individual error bars. In SA 93, a proper 
motion offset was apparent between faint (g > 19.0) candidates and brighter 
Sgr candidates; for this reason, only the faint candidates were retained. The 
red asterisks represent the final maximum likelihood estimate of the mean 
Sgr debiis motion in each field. 

Sgr MSTO candidates for follow-up spectroscopy. Because 
the available proper-motion data do not sample the MSTO as 
robustly as in SA 94, this field was given lower priority for 
spectroscopy, with only one relatively short MMT-i-Hectospec 
configuration observed (see Table |2]l. However, in this one 
Hectospec setup, nearly all available g < 20.2 Sgr candidates 
satisfying the (somewhat relaxed relative to S A 94) target 
selecti on criteria were observed. The overlaid [Girardi et al.l 
(12004ft isochrone for an old (10 Gyr), metal-poor ([Fe/H] = 
-1.3) population at the expected distance {d = 28 kpc) of Sgr 
debris in SA 93 (blue curve in Figure [TTli is seemingly con- 
sistent with all but the brightest of the identified candidates. 
We have also overplotted (dashed gray curve) an older (^15.8 
Gyr), more metal-poor ([Fe/H] = -1.7) isochrone with the 
same distance as the blue curve; this curve passes through the 
positions of the bright, blue stars, suggesting that they may 
also be Sgr members of an old horizontal-branch/asymptotic 
giant branch population. However, the MSTO of such an 
old population is beyond the magnitude limit of our proper- 
motion survevFl The RPMD (Figure [T3T l for SA 93 shows 
the same ridgeline as in the CMD, shifted by the final mea- 
sured tangential velocity for the Sgr candidates. Because the 
uncertainty in the mean proper motion is much larger in this 
field than in SA 94 (on the order of '--^lOO km s~' in tangen- 
tial velocity for SA 93, compared to ^-^50 km s"' in SA 94), 
it is difficult to conclude much from the RPMD. The benefits 
of the additional 4-meter plates in SAs 94 and 71 are clearly 
illustrated by the relatively fewer identified Sgr stream mem- 
bers in the shallower SA 93 field compared to the deeper data 
sets. Zoomed-in versions of the VPDs for SAs 94 and 93 are 
given in Figure [14] with error bars shown on all points to il- 
lustrate the quality of the proper-motion data (note that the 
error bars on individual stars in each of the panels of Fig. [14] 
are of comparable size, in spite of the higher-quality proper 
motion data in SA 94 than in SA 93, because the majority of 
Sgr candidates in SA 94 are faint [g > 21] MSTO candidates, 
while SA 93 candidates are mostly > 1 magnitude brighter 
than this). As previously mentioned, a proper motion offset 
was seen between faint {g > 19.0) stars in SA 93 and brighter 
candidates; we included only the faint (g > 19.0) stars in our 
proper motion measurement, because these are more likely to 
be true Sgr members. The maximum likelihood estimate of 

Note that we are not suggesting that Sgr has stars 15.8 Gyr old - this 
isochrone is simply meant as a guide to show that these stars could plausibly 
be BHB stars associated with Sgr. A blue horizontal branch is only seen 
in the oldest (log(age) > 10.15) of the iGirardi et al.ll2004l isochrones at this 
metallicity. 



the absolute proper motion of Sgr debris in each field is rep- 
resented by the large red asterisks, which have la uncertainty 
smaller than the size of the point. 

3.2.3. Summary of Sgr Candidate Selection 

In summary, the basic Sgr candidate selection in each Se- 
lected Area began with a broad RV selection centered on the 
apparent Sgr velocity peak (e.g., -250 < Vhei < -100 km s"' in 
SAs 94 and 93). This was followed by removing high proper 
motion stars (|/i| > 10 mas yr"' in either dimension), which 
would have tangential velocities much greater than the Milky 
Way escape velocity if those stars were at the 25-40 kpc dis- 
tances of Sgr debris along the SA lines of sight. We then used 
our knowledge of the distance and metallicity expected for 
Sgr debris in these fields to remove stars at faint magnitudes 
that are more than ^0.3 magnitudes from the expected color- 
magnitude and RPMD loci in each field (note that in some 
cases we did not apply this at bright magnitudes, since there is 
considerable uncertainty about the exact CMD locus for Sgr 
debris above the turnoff). We calculated mean proper mo- 
tions from the selected stars, then iteratively removed outliers 
(> 3a) in proper motion. Finally, we manually inspected the 
remaining candidates to remove any stars that were ~ 1 - 2cr 
outliers from the identified Sgr locus in all observables (i.e., 
RV, g-r color, magnitude, and proper motion). 

Once Sgr candidates were selected using color, magnitude, 
RV, proper motion, and RPMD criteria, the kinematical prop- 
erties of Sgr debris in each SA field were estimated using 
a maximum likelihood method (e.g., jPryor & Mevlan 1993c 
iHargreaves eTal] 1 1 994 iKlevna et alj|2002h . The final mea- 
sured radial velocities and velocity dispersions, along with 
uncertainties in these values, are given in Table [3] in both the 
heliocentric and Galactocentric (GSR) frames. The estimates 
of the mean radial velocity (in the GSR frame) of Sgr de- 
bris can be seen in the top panels of Figure [15] which over- 
lays our measurements a top the best-fitting debris model of 
iLaw & Maiewskil ( l2010al) . T he colors in the figure wer e cho- 
sen to match those used by iLaw & Majewskil (l2010al) . who 
color-coded points along the stream according to the orbital 
passage on which they became unbound. Gold points corre- 
spond to debris stripped during the two most recent perigalac- 
tic passages, and magenta during the previous two passages. 
This portion of the trailing tail has a well-constrained radia l 
velocity tren d, which was measured bv lMaiewski et alj (120041) 
and Monaco et al. (2007), and was one of the constraints on 
the Law & Majewski (2010a) models. Sgr model trailing tail 
debris within ±3.0° of the position of each Kapteyn line of 
sight is shown in the right-hand panels of Figure [15] as small 
open squares. This illustrates that not only do our measured 
radial velocities agree quite well with the model, but that the 
dispersion in each field appears to match the distribution of 
model points. However, we caution that the velocity disper- 
sions we find (see Tab le [3] ) are higher than thos e derived by 
iMaiewski et al.l (120041) and lMonaco et alj (l2007h for the trail- 
ing tail. It is unclear whether the dispersion is truly intrinsi- 
cally higher for Sgr trailing tail MSTO populations we sam- 
pled than for the M giants previously studied, or whether our 
derived dispersions are inflated by the large measurement un- 
certainties for our RVs, velocity zero-point offsets between 
the many data sets we have combined, or Milky Way fore- 
ground/background contamination in our Sgr candidate sam- 
ples. While we have endeavored to account for all of these 
factors, a robust conclusion likely requires high-resolution 
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spectra of trailing-tail MSTO stars. 

3.3. Proper Motions 

With the sample of Sgr debris candidates identified in each 
field, we used a maximum-likelihood method to estimate 
the Sgr absolute proper motions in both spatial directions 
(HaCosiS), fis). These results are seen in Table [3] with their 
uncertainties, which include the uncertainty in the proper mo- 
tion zero point in each field added in quadrature to the max- 
imum likelihood error estimate. The uncertainties for each 
field depend on many factors, including the depth and quality 
of the plates, the number of available reference objects (i.e., 
background QSOs and galaxies) used to convert from relative 
to absolute proper motions, the number and depth of spec- 
troscopic targets obtained (and thus the number of Sgr candi- 
dates identified), and the position of each field relative to the 
highest-density regions of the stream (i.e., the number of Sgr 
candidates expected in each field). 

Mean proper motions along Galactic coordinates (i.e., 
UiCOsb,iJ,h) in each fiel d are compared to the model of 
iLaw & Majewskil (|20 10a') as a function of A© in the middle 
and lower panels of Figure [15] with points once again color- 
coded by the orbital passage in which they became unbound, 
and model points within ±3.0° of each SA line of sight high- 
lighted for guidance. The results for SAs 94, 93, and 117 
agree nicely with the model predictions within the uncertain- 
ties in III,, but show a ^ I -2a offset from the main trend 
in fiicos{b) (left middle panel). The mean proper motions 
along both directions for SA 71 (the leftmost data points in 
Figure [ISll are slightly shifted (by '--^ 1 - 1. 5a) from the mean 
of the model prediction for this field. A number of factors 
contribute to the difficulty in selecting a "pure" sample of Sgr 
debris in SA 71, and there is thus an additional uncertainty 
(besides the formal errors) in the mean Sgr debris proper mo- 
tions in this field. First, this field is at somewhat low lati- 
tude (b = -34.7°), and thus suffers greater contamination from 
Galactic populations, as evidenced by the extended tail of 
low-velocity stars in the lower panel of Figure |5] This lower 
latitude also means that S A 7 1 suffers significantly more red- 
dening than higher-latitude fields - nearly 0.6 magnitudes of 
extinction using the Schlegel et al. (1998) maps. Secondly, 
the distance of Sgr debris increases with A© along the por- 
tion of the trailing tail in this study, so that stream members 
in SA 7 1 are nearly 40 kpc away, making them fainter than in 
the other fields. Finally, the poor-quality photographic pho- 
tometry we are limited to in this field renders inscrutable the 
typical CMD features such as the blue edge of the Galactic 
disk MSTO and the Sgr upper main sequence. 

Assuming distances to Sgr debris in each field as given in 
Tableland values of Vciic = 220 km s"' and Rq = 8.0 kpc, 
we converted the measured Sgr debris motions to Galacto- 
centric UVWcc velocities (i.e., Cartesian velocities such that 
the Sun is moving at 9.0, 232.0, 7.0 km s~' assuming Vcm 
= 220 km s"' and (Up, Vp, Wo) = (9.0, 12.0, 7.0) km s"' 
jMihalas & Binnevlll981b for the solar motion relative to the 
Local Standard of Rest; to facilitate direct comparison, values 
used for these constants are the same as those in the ILMIOI 
model), which are shown in Table|3](note that we placed SAs 
92 and 116 - the two fields with no securely identified Sgr 
members - in a separate section in Table [3] These data are 
given for completeness, but are not used for subsequent anal- 
ysis.). The Ucc component should dominate the total space 
velocity of Sgr debris in each of these trailing arm fields. This 



can be surmised from Figure |2] which shows that the mo- 
tion of the Sgr trailing tail is oriented almost parallel to the 
Galactic X-axis in the Xqc-Zgc plane. This, in addition to 
the fact that the Sagittarius orbital plane is only slightly mis- 
aligned with the Galactic Xgc-Zqc plane (Mai ewski et al] 
2003, 2006), suggests that most of the motion in this part 
of the trailing tail is inward toward the Galactic center and 
roughly parallel to the Galactic disk (i.e., the Xgc-Jgc plane) 
at a distance of ^^20-25 kpc below the disk. This is borne 
out by the measured UVWgc velocities in our SA fields - the 
Ucc component in the four fields with quality measurements 
is by far the largest component of the 3-D motion. This can 
be seen even more clearly by considering the proper motions 
along Galactic coordinates, but in a Galactic rest frame (des- 
ignated as /i{cos(b) and /i^). These proper motions, given 
in Table [H show cos(b) proper motions of nearly zero in 
each field - as expected for streaming motions confined to the 
Xgc --Zgc -plane. As we will show in Section 4, the offset 
of these derived longitudinal proper motions (reflected in the 
Vgc Galactic velocity component) from zero can be used to 
reevaluate the velocity of the Local Standard of Rest (under 
the assumption that the longitudinal motions should be zero). 

4. CONSTRAINTS ON MILKY WAY STUCTURE 

As discussed in Section 1.1, the opportune orientation of 
the Sgr trailing tidal tail means that the observed motion of 
tidal stream stars in the Galactic Y direction (i.e., towards 
{l,b~\ = [90°, 0°]) is dominated by the solar reflex motion, 
which consists of the solar peculiar motion and the Galactic 
rotational motion at the solar c ircle ( i.e. , the Local Standard 
of Rest 8lsr)- As shown by Maiewski et al] ( |2(K)6|) . the in- 
trinsic motion of Sgr debris along the Y direction (Vgc, con- 
tained primarily in the /// cos(b) component of proper motion) 
varies only slowly across the region of the trailing tail be- 
tween 70° < Aq < 130°, making the fields of view in which 
we have deep proper motion data ideal for constraining 8lsr- 
It can be seen in Table [3] that Vgc for Sgr debris in each of 
the four fields (SAs 71, 94, 93, and 1 17) with reliable data is 
non-zero at the ^ \a level. Setting aside SA 71, in which it 
is difficult to securely identify Sgr debris, the remaining three 
fields exhibit Vgc systematically offset to negative values. If 
indeed the expected Vgc for Sgr debris in these fields is zero, 
this suggests that the value of 0lsr that was subtracted from 
the V-component of these velocities was lower than it should 
be - i.e., Glsr should be greater than the canonical 220 km 
s-i. 

In this section, we use variations on the iLMld numerical 
model of the Sgr tidal stream to isolate the contribution to 
jii cos(b) from the solar reflex motion and identify the value 
of Olsr favored by our proper motion data. 

4.1. N -body models 

Though the measured Galactic Cartesian V-velocity (i.e., 
motion along the Galactic F-component) of Sgr trailing tidal 
debris is dominated by Solar reflex motion, there is some con- 
tribution of intrinsic Sgr motion to the V-component of debris 
velocities. In particular, we must consider the following ef- 
fects when trying to back out Glsr from measured Vgc for 
Sgr debris: (1) the slight inclination of the Sgr debris plane to 
the Galactic XZgc plane means a small fraction of Sgr space 
motion is projected onto the measured motions (i.e., the V 
velocity is in fact a function of both 8lsr and intrinsic Sgr 
motion); (2) the Galactic Standard of Rest (GSR) frame ra- 
dial velocities used to constrain the Sgr model were derived 
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Table 3 

Sagittarius Stream Kinematics in Kapteyn Selected Areas 



SA 


N 


^helio 




O"0 


fla COs(5) 




Hi cos(b) 




Uac 


Vac 


Wgc 


distance 






km s"' 


km s"' 


km s"' 


mas yr"' 


mas yr"' 


mas yr"' 


mas yr"' 


km s"' 


km s"' 


km s"' 


kpc 


71 


33 


-172.9±4.9 


-141.4±4.9 


25.6±3.5 


0.39±0.44 


-0.55±0.42 


0.66±0.43 


-0.17±0.43 


138.1±47.1 


80.7± 77.1 


80.1±63.3 


38.0±4.5 


94 


64 


-142.4±2.3 


-141.3±2.3 


17.4±1.8 


0.35±0.27 


-2.32±0.29 


1.88±0.28 


-1.40±0.28 


227.9± 34.1 


-49.5± 54.0 


-12.6± 30.8 


29.5±4.0 


93 


15 


-157.2±4.7 


-114.1±4.7 


16.1±3.3 


0.42±0.73 


-2.66±0.56 


1.60±0.70 


-2.16±0.60 


210.6± 81.3 


-101.6± 88.4 


-10.7±45.4 


28.0±3.0 


117 


43 


-90.0±3.9 


-69.3±3.9 


24.0±2.8 


0.02±0.64 


-3.11±0.72 


1.28±0.65 


-2.84±0.71 


229.3± 83.7 


-76.8± 94.6 


11.5±24.9 


25.0±4.0 


92 


6 


-158.9±6.3 


-78.4±6.3 


10.6±8.7 


1.34±1.07 


-3.84±0.66 


1.46±1.()6 


-3.80±0.66 


145.6±123.3 


-296.9± 107.9 


-84.6± 45.6 


27.5±2.5 


116 


10 


-75.4±6.7 


-22.3±6.7 


20.3±4.9 


1.12±1.04 


-3.30±1.15 


-0.67±1.07 


-3.42±1.12 


87.3±124.3 


-170.8±129.6 


-22.7± 35.3 


24.5±2.0 



Note. — All calculations assume Vciic 
(in a right-handed frame). 



= 220 km s"' at Rq = 8.0 kpc. We used the solar peculiar motion of lMihalas & BinneM fT98Ti) : (Uo, Vq, Wq) = (9.0, 12.0, 7.0) km s" 
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Figure 15. Kinematics of Sgr candidate stars measured in four SA fields (SAs 71, 94, 93, and 117, from left to right) as a function of longitude in the Sgr 
coordinate system defined by Maiewski et al. (2003). From top to bottom, the panels depict GSR-frame radial velocity (Vgsr), proper motion along Galactic 
longitude (pi cos b), and proper motion along latitude (pi,)- Colored points depict trailing tail debris from the best-fit Sagittarius model of Law & Majews^ 
(2010a), with different colors representing debris stripped on successive orbits, as in lLaw & Majewski (2010a). In each row, the right-hand panels depict stars 
selected from the model to be within ±3° of each of the SA fields as small black squares; these are the points to which we compare the measured kinematics 
in each field. In the left panels, large open diamonds with error bars represent the maximum likelihood estimates of the mean kinematics of Sgr debris. Note 
that it is important to compare measured kinematics (on the left) to the "clouds" of small black squares in corresponding fields in the right panels, rather than 
comparing to the trend defined by all of the colored points. This is necessary because at a given Aq, the model debris plotted here can come from a large area on 
the sky. We are interested in comparing only what the model predicts "should be" seen in each pencil-beam field of view, and thus we select only model debris in 
corresponding regions of the sky. The measured Vqsr for Sgr debris matches the model very well in all four fields, and the fij, proper motions agree in three of 
the four fields (see the text for discussion of the difficulties in selecting Sgr debris in SA 71, the ~ Icr discrepant point furthest to the left). In /i; cos b, however, 
the measurements for all but one field (again, SA 71) are systematically offset to higher proper motions than predicted by the model. We show that this offset can 
be accounted for by an upward revision of 0lsr. ths rotation speed at the Solar circle. 



assuming a value of 6lsr; (3) changing 6lsr correspond- 
ingly changes the Milky Way mass scale, which thus affects 
the space velocity of the Sgr dSph in the models. Therefore, 
taking into account th ese dep endencies on the assumed value 
of Qlsr, we repeat the LMIO analysis, changing 6lsr to con- 
struct self-consistent models for the Sgr tidal stream in each 
of four choices for the Local Standard of Rest speed, namely 
©LSR = 190,250,280,310 km s"' (in addition to the original 
|LM10 value of 8lsr = 220 km s"' used in earlier sections of 
this paper). 



Our methodology is described in detail bv lLMlOi In brief, 
we constrain the model Sgr dwarf to lie at the ob served loca- 
tion (^ / ^)= (5.6°, -14.2° ), distance Dsg, = 28 kpc (I Siegel et al.1 
1201 ll ISiegel et al.ll2007b . and radial velocity vsgi = 142.1 km 
s"' in the heliocentric frame. The orbital plane is constrained 
to be that defined by the trailing arm tidal debris, which 
has e xperienced minimal angular precession ([Johnston et al.l 
l2005h . and the speed of Sgr tangential to the line of sight (vtan; 
see Section 3.3 of LMIO) is constrained by a minimization 
fit to the radial velocities of trailing arm tidal debris. We fix 
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Table 4 

Galactic Frame-of-Rest Proper Motions of Sagittarius Debris in Selected 
Areas 



SA 




fj,'^ cos((5) 


^5 


fi'l cos(b) 






(kms-') 


(mas yr"') 


(mas yr ' ) 


(mas yr"' ) 


(mas yr ' ) 


71 


-141.4 


-0.47 


0.39 


-0.61 


0.00 


94 


-141.3 


-0.77 


-1.09 


0.22 


-1.31 


93 


-114.1 


-0.71 


-1.36 


0.00 


-1.54 


117 


-69.3 


-1.17 


-1.57 


-0.44 


-1.91 



"All calculations assume V^irc = 220 km s ' at = 8.0 kpc, and solar pe- 
culiar motion of {Uo, Vq, Wq) = (9.0, 12.0, 7.0) km s"' (in a right-handed 
frame). 



Table 5 

Masses of Galactic Bulge and Disk Components in the Sagittarius Models 



©LSR 






a 


(kms-') 


(Me) 


(Mq) 




190 


6.8 X lO'O 


2.3 X lO'" 


0.68 


220 


1.0 X 10" 


3.4 X 10'° 


1.00 


250 


1.4 X 10" 


4.6 X 10'° 


1.35 


280 


1.8 X 10" 


6.0 X 10'° 


1.76 


310 


2.2 X 10" 


7.5 X lO'O 


2.21 


264 


1.5 X 10" 


5.2 X 10'" 


1.53 


232 


1.1 X 10" 


3.9 X 10'° 


1.14 



"The mass of the disk and bulge components in each of the models of the 
Sgr stream. Each model is specified by the value of 0lsr that constrained 
the fit; die 220 km s"' model is that of lLMlfl The ratio of disk to bulge 
mass is constant throughout - the constant a is the scaling factor, such that 
Mdisk = a * Mjisk.220km/s and Mb„ige = a* M|,„,ge,220km/s ■ The total mass, axis 
ratios, and scale length of the Galactic dark matter halo were fixed to the 
best-fit values of lLMlOt 

the mass and radial scalelength of the Sgr progenitor so that 
the fractional mass loss history of the dwarf is similar in all 
models to that of LMIO. 

The adopted Milky Way Galactic mass model consists of 
three component s: a Hernquist spheroid ( representing the 
Galactic bulge), a lMivamoto & Nagail ( 119751) disk, and a log- 
arithmic dark matter halo. The Local Standard of Rest in this 
model is given by: 

©LSR = ^/^©(abulge + adisk + ahalo) (1) 

where abuige. Odisk. and Ohaio respectively represent the grav- 
itational acceleration exerted on a unit-mass at the location 
of the Sun due to the Galactic bulge, disk, and halo compo- 
nents. In the ILMIOI model (for which Glsr = 220 km s"'), 
the bulge/disk/halo respectively contribute 32%/49%/19% of 
the total centripetal acceleration at the position of the Sun, 
corresponding to bulge/disk masses Mbuige = 3.4 X lO'^Mo 
and Mdisk = 1.0 x 10"Mo, and a total mass within 50 kpc of 
4.5 X 10"Mo. 

Since the baryonic Galactic disk and bulge components are 
the dominant factors in determining 6lsr (together compris- 
ing > 80% of the total centripetal force), we therefore scale 
the total bulge -I- disk mass as necessary to normalize the rota- 
tion curve at the solar circle (Rq = 8 kpc) to the chosen value 
of Qlsr- The masses of the disk and bulge components in 
each of the models are given in Table |5] We leave the ra- 



tio of bulge/disk mass fixed in order to preserve the shape of 
the rotation curve interior to the solar circle. In addition, we 
fix the Galactic dark matter halo parameters (mass, axis ra- 
tios, and scalelength) to the best-fit values derived by LMlO 
since these authors found that these values were relatively in- 
sensitive to factors of 2 variation in the mass scale of the 
baryonic Galactic components. 

4.2. Results 

Constraints on 9lsr were derived in two ways. In the first 
of these methods, we assumed (as argued previously in this 
paper, as well as in MLPP) that the dominant contribution to 
the measured /i/ cos(b) component of Sgr motion is due to the 
solar rotation, and that the largest component of Glsr is along 
/i; cos(b). We have shown that these are reasonable first-order 
assumptions, and thus use only the longitudinal proper mo- 
tions as constraints on fitting 8lsr in our first attempt. After 
doing so, however, we performed a similar analysis, but us- 
ing all three dimensions of Sgr debris motions as constraints 
to determine 0lsr- In the following, we present both results, 
which come out somewhat different from each other (though 
consistent within la). Fits using only /i/cos(b) tend to prefer 
relatively high values of 0lsr, while those constrained by full 
3-D kinematics tend toward lower values more in line with the 
lAU standard of 220 km s"' . 

4.2.1. Qlsr Constraints Using Only /i;COs(Z?) Motions of Sgr 

Debris 

We quantify the agreement of our proper motions with 
those of simulated Sgr tidal debris from each of our grid of 
models using a statistic. The fitting was performed 
using mean Sgr debris proper motio ns in only SAs 71, 94, 
93, and 117 - as discussed in Section [TTl the results in SAs 
92 and 116 are unreliable for a vari ety of reasons. For our 
model comparison, we first select all lLMlOl Sgr model points 
within ±3.0° in both Aq and (a,S) of each SA field. The 
large area (relative to the 40' x 40' coverage of each SA field) 
used to select model debris corresponding to each S A position 
ensures that enough A^-body particles are selected for robust 
measurement of model debris motions at each position. This 
also makes the fitting less sensitive to small-scale differences 
in positions and densities of debris stars between the mod- 
els and the actual stream that arise due to the vagaries of the 
modeling and our incomplete knowledge of the Sgr trailing 
tail properties. Figure [16] shows the model debris /i/cos(b) as 
a function of A0 for each of the five Sgr simulations, with 
points corresponding to each SA field shown as small open 
gray squares. It is clear that /i/cos(b) changes very little over 
the 6° ranges in (a, S) used. Furthermore, the small number of 
selected model points, even in such a large selection region, 
shows that these broad selection criteria are necessary to have 
sufficient model points for comparison. The maximum like- 
lihood proper motion results in SAs 71, 94, 93, and 117 are 
shown in Figure [16] as open black diamonds, with error bars 
reflecting Icr uncertainties. It can be seen in the figure that the 
models with 8lsr > 220 km s"' tend to reproduce the longi- 
tudinal proper motions for most of the fields better than the 
standard 220 km s"' value of this fundamental constant. Note 
that SA 71, at Aq = 128°, is the exception to this trend - as 
discussed in Section 3.3, identifying bona fide Sgr debris in 
this field is more difficult than the others and so these data are 
more suspect. 

Each of the Sgr model simulations provides predicted kine- 
matics of Sgr debris for a Galactic potential constrained by a 
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Figure 16. Mean longitudinal proper motions, /^;cos(b) (large open diamonds with eiTor bars), in (from left to light) SAs 71, 94, 93, and 117 as a function of 
Aq. Model debris from the Sgr simulations are shown as the colored points; from top to bottom, these represent models with ©lsr = 310, 280, 250, 220, and 
190 km s"' . The sudden drop in fii cos(b) for Aq < 90° is due to the inversion in sign that occurs as the debiis sweeps past the South Galactic pole. Small open 
gray squares denote the model debris corresponding (within ±3° in RA, Dec, and Aq) to each SA field. It is clear that higher values of 0lsr provide a better 
match of the small grey squares to the observed /j; cos(b) values in the Selected Areas (diamonds). 



given value of 8lsr- As discussed previously, the Galactic V- 
component of the Sgr debris space velocity along the trailing 
tail contains little contribution due to the Sgr motion; nearly 
all the V velocity (as measured by the /i/ cos(b) component 
of the proper motion) is reflected Solar motion. Thus, to first 
order, we can simply compare our mean Sgr proper motions 
along Galactic longitude in the trailing tidal tail to our models 
of Sgr debris for different values of Glsr and determine the 
value of 0LSR that best reproduces the measured PMs. To do 
so, we defined a residual: 

2 _ M/.SA[i]-Ail,mod[i] 



where /i/.sA[i] represents the mean /i/cos(b) proper motion 
in each of the four SA fields, and /x/ niod[i] the mean proper 
motion of the corresponding model debris for each field. 



The residuals are weighted by the uncertainty, cr^,.sA[i], in 
each SA proper motion. This statistic was initially cal- 
culated for the proper motions presented in Table [3] rela- 
tive to each of the five Sgr debris models (corresponding to 
6lsr = 190,220,250,280,310 km s"'). Rather than running 
new (and laborious) A^-body simulations for many interme- 
diate values of 8lsr and calculating x^ for each of them, 
we choose to find the minimum x^ by fitting a parabola to 
the x^ results for each of the five modeled values of 8lsr- 
The results of the calculation and the parabolic fit are seen 
in Figure [17] for the proper motions of the "best" SA sam- 
ples given in Table [3] The minimum of the parabola yields 
BLSR.min = 270.4 km s"'. 

To estimate the uncertainty in 8lsr, we c hoose a boo tstrap 
(resampling with replacement) method (see I Andra3l20 1 and 
references therein). This technique uses the entire sample of 
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individual Sgr candidate star proper motions, and thus yields 
an estimate of the errors in 8lsr including the effects of 
proper motion measurement errors and "contamination" of 
the proper-motion samples by Milky Way stars. Our selected 
samples of Sgr candidates in SAs 71, 94, 93, and 117 should 
contain mostly Sgr debris, plus some amount of contamina- 
tion by MW stars that will vary depending on the depth of 
the proper-motion and radial-velocity catalogs, the local Sgr 
stream density, and the Galactic latitude of each field. From 
the Sgr samples in each field, we performed 100,000 boot- 
strap resamplings, wherein random selections were made 
from the original stars in each field (i.e., the catalogs of 
candidates were resampled with replacement). Iteratively 3(t- 
clipped mean proper motions of these resampled Sgr candi- 
dates were measured, and the mean proper motions used in an 
identical fitting routine to that described above. Assum- 
ing that the contaminants in each sample are somewhat uni- 
formly distributed in their kinematical quantities, this method 
should yield a statistically robust result for 0lsr and its un- 
certainty (due to both the intrinsic measurement errors and the 
MW contamination). The best-fitting values of 9lsr for these 
100,000 samples are given as a histogram in the left panel of 
Figure[T8] along with a Gaussian fit (red curve) to the results. 
From this Gaussian, we derive a final value of 8lsr = 264 ±23 
km s"'. 

4.2.2. 0LSR Constraints Using Three-Dimensional Motions of Sgr 

Debris 

The constraints we derived on 0lsr using only the longi- 
tudinal proper motions assume that the contributions of 0lsr 
to Hb and Vgsr are negligible. If the Sgr orbital plane was 
exactly coincident with the Galactic XZqc plane, then indeed 
the rotation velocity at the solar circle would only be reflected 
in the /i/cos(b) motions. In reality, the Sgr orbital plane is not 
perfectly aligned with the Milky Way XZ-plane, so there is 
some projection of Vdrc onto ^i, and Vgsr- In fact, for Sgr de- 
bris in the four fields of view comprising this study (SAs 71, 
94, 93, and 1 17), only (75%, 87%, 61%, 58%, respectively) of 
the total value of 6lsr is projected onto cos(b). We ran the 
fitting again, but this time including all three dimensions 
of the motion as constraints. The bootstrap analysis gave a re- 
sult of 0LSR = 232 ± 14 km s"' - a histogram of the bootstrap 
results is seen in the right panel of Figure[T8] This mean value 
is lower by ^ 1.4cr than the result using only /i/cos(b). For- 
mally, this is a better fit than the one-dimensional result (with 
uncertainty of 14 km s"' compared to an uncertainty of 23 km 
s"' from the fits using only longitudinal proper motions), but 
the two are consistent within their la uncertainties. 

Finally, we performed the same exercise using all three di- 
mensions of Sgr debris motions, but excluding the less reli- 
able SA 71 field. The uncertain identification of Sgr debris in 
SA 71 is likely the reason this field (at A© = 128°) is an out- 
lier from the predicted kinematical trends in Figures [15] and 
[T6l The bootstrap fit using only SAs 94, 93, and 117 yields 
©LSR = 244 ±17 km s"' . This slightly higher value for Glsr 
suggests that (as is evident in Fi gures [TS] and [T6]l the Sgr can- 
didates in SA 71 skew our results toward lower Glsr- 

Ultimately, we have derived three estimates of 6lsr - one 
based on a simple one-dimensional analysis (using all four 
fields) that gave 264 ±23 km s"', another based on three- 
dimensional data yielding 232 ± 14 km s"', and a final 3-D 
result with SA 71 excluded, which gave 244 ± 17 km s"'. It 
is likely that the true result is somewhere between the two 
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Figure 17. Total residuals for the final mean proper motions relative to 
corresponding model debris. Each point represents a for one of the five 
models in which we vary ©lsr from 190-310 km s"'. A parabola fit to the 
results (the red curve) yields a minimum x^ at 0lsr = 270.4 km s"'. 





Figure 18. Resulting ©lsr corresponding to the minimum x of 100,000 
bootstrap resamplings of the individual stellar proper motions in each SA 
field. The left panel shows fits using only the /i/ cos(b) component of Sgr 
debris motions in SAs 71, 94, 93, and 1 17. A Gaussian fit (overlaid as the red 
curve) to the results yields ©lsr = 264 ± 23 km s"'. The right panel shows 
the results using all three dimensions of Sgr debris motions (i.e., /i;cos(b), 
fil,, and Vgsr) ^s constraints, which yields ©lsr = 232 ± 14 km s"' . 

extremes (264 km s"' and 232 km s"^) from our methods. 

4.2.3. Sgr Disruption Models For Best-Fitting Slsr 

We now repeat the A^-body analysis described in Section|4T| 
two times, first taking 8lsr = 264 km s~' as a constraint on 
the models, then again using Olsr = 232 km s~'. The result- 
ing A^-body model for the 264 km s"' case matches the angular 
position, distance, and radial velocity trends of the observed 
Sgr tidal streams (using all of the observational constraints in- 
cluded in the original LMIO model) almost equally as well as 
did the LMIO model (formally, x = 3. 1 for the 6lsr = 264 km 
s"' model, compared to x = 3.4 for the LMIO model; see dis- 
cussion in Section 4.3 of LMIO). In addition, as demonstrated 
in the left panel of Figure[T9]the proper motion of the remnant 
core of the Sgr dwarf in this revised model (/i/cos(b)= -2.54 
mas yr"', /i^ = 1.92 mas yr~') is a s ubstantially better match 
to observations (e.g., Din escu et al . 2005; Prvoretal. 201(J 
than was the LMIO model. 

However, the A^-body model in a Milky Way halo with 
Blsr = 232 km s"' fits equally well as does the 264 km s"' 
case, X = 3 . 1 . For this model, the Sgr core proper motion (seen 
in the right panel of Figure [T9]l is intermediate between those 
of the LMIO model and the 264 km s"' result, as might be 
expecte d. In this case, the pro per motions are d i screpa nt with 
both the lDinescu et all ( l2005h and iPrvor etall ( 120 1(3) results 
at the ~ 1 .5(7 level. In the following subsection we discuss the 
ramifications of these 0lsr results for the Milky Way halo. 
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Figure 19. Proper motion estimates for tlie Sgr core in Galactic coordinates. 
Blue and red shaded ellipses show 1 u and 2(t uncertainty regions around the 
measurements of Dinescu et al. (2005) and Pryor et al. (2010) respectively. 
The proper motion of the Sgr model dwarf described by LMIO (©lsr = 220 
km s"') is indicated by a filled black square. In the left panel the Icr range 
of proper motions corresponding to the value ©lsr = 264 it 23 km s"' we 
found using only the /.t; cos(b) component of Sgr debris motions is indicated 
by the eiTor bars surrounding the filled black circle. The right panel is similar, 
but for the result (©lsr = 232 ±14 km s"' ) using all three dimensions of Sgr 
debris kinematics. The orientation of the error bars represents the direction to 
which changes in ©lsr correspond in this diagram; the proper motion of the 
Sgr core is constrained to lie along this line because the orbital plane is fixed 
by the observed position of the tidal debris leading and trailing Sgr, which 
trace its orbit. 

4.3. Robustness of the Galactic Mass Models 

The results for 0lsr from our analysis in the previous sub- 
section suggest that the true value of 8lsr as constrained by 
Sgr trailing tail debris likely lies between 232-264 km s~' , but 
that more (or more sensitive) proper motion measurements of 
Sgr trailing debris are needed to resolve this issue. In this 
subsection, we will discuss the implications of the 6lsr con- 
straints resulting from our two methods; we remind the reader 
that the upper end of the range (i.e., the 264 km s"' result) is 
less robustly determined than the results that produced lower 
values of the circular velocity. However, we include this value 
in our discussion to present the reader with the range of possi- 
ble ramifications of what, in either case, represents an upward 
revision of 6lsr from the accepted value. 

The value of &lsr = 232 ± 14 km s"' we found using all 
three dimensions of Sgr debris kinematics is consistent with 
the canonical 220 km s"' value at the roughly la level. How- 
ever, this value is also consistent (within the Icr uncertain- 
ties) with recent determinations of 6lsr that have found the 
rotation spe ed to be higher than the standard 220 km s"' 
value [e g., IReid et aH (i009h - Olsr = (254 ± 16)(7?o/8.4) 
km s-i ; iBovv etal.l (l2009h - Glsr = 244 ± 13 km s"!]. For 
a change in 9lsr of only about 10 km s"', it is difficult to 
make any conclusions about whether the additional mass re- 
quired to increase the rotation speed must reside in the Galac- 
tic halo or the disk/bulge. We note that placing the addi- 
tional mass in the disk and bulge (with the halo fixed) yields 
Mbuige = 3.9 X 10'°Mo and Mdisk = 1.1 x 10 "Mq for the 232 
km s~' model - an incre ase of ^ 10% over the disk and bulge 
mass from the model of iLMlOt The relatively high value of 
©LSR = 264 ± 23 km s"' found by our analysis using only the 
/i/ cos(b) motions would require that the mass of the Galac- 
tic bulge and disk components be increased by ~ 50% from 
the values assumed bv lLMlOl to Mbuige = 5.2 x lO'^M© and 
Mdisk = 1.53 X lO'^M©. This disk mass is near the peak of 
the probability distribution (Mdisk = 1 .35 x 1O"M0) found by 
iKoposov et al.l (^010) based on fitting the GD-1 stream in a 
three-component gravitational potential similar to our own. 
We caution however that the orbit of Sgr is largely insensi- 
tive to the distribution of the excess mass between the two 
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Figure 20. Model Milky Way rotation curves as a function of radius from 
the Galactic center Solid blue/green/red lines respectively represent rotation 
curves with ©lsr = 264 km s"' achieved via scaling the Galactic bulge+disk. 
Galactic disk alone, and Galactic halo alone. Included for comparison is the 
rotation curve of the original LMIO model (solid black line) normalized to 
©LSR = 220 km s"' . The vertical dotted fine represents the location of the 
Sun at Rq = 8 kpc. 

baryonic components, and solutions that yield similar can 
be found by ascribing all or part of the needed adjustment in 
8lsr to changes in the mass of either the disk or bulge com- 
ponents alone. We do note, however, that the relative fraction 
of the total disk-i-bulge mass in each component is constrained 
by the need to reproduce the shape of the observed Milky Way 
rotation curve interior to Rq (see Figurel20]l. 

The total mass of the Milky Way interior to 50 kpc in 
the 264 km s"' model is 5.2 x 10"Mq, similar to the value 
of 4.5 X lO^M© in the iLMld model. Since we have ac- 
counted for the increased 8lsr = 264 km s~' by increasing 
the disk-i-bulge mass (which is a relatively small component 
of the total virial mass), the mass of the Milky Way interior 
to 200 kpc is Mvii- =1.6 x IO'^M q, similar to the value of 
1.5 X 10'^Mq derived by ILMIOI assuming that Olsr = 220 
km s"^ . 

We note that it was not possible to obtain a satisfactory 
model (within our parameterization of the Milky Way com- 
ponents; exploration of different dark halo models is beyond 
the scope of this work) for the Sgr stream by leaving both the 
bulge and disk masses fixed at their iLM lOl value s and account- 
ing for changes in 8lsr by scaling the dark matter halo. The 
dark matter halo profile is characterized by the parameters 
I'haio and rjiaio (see Eqn. 3 of LMIOQ which describe the total 
mass normalization and radial scalelength of the halo respec- 
tively. Since dark matter in the LMIO model contributes only 
19% of the total centripetal acceleration in the solar neighbor- 
hood, Vhaio must be scaled up drastically (by a factor of ^ 3 in 
total halo mass) to increase Olsr from 220 km s"' to 264 km 
s"', and necessitates a large increase in the space velocity of 
Sgr along its orbit (to ^ 400 km s"') to produce a leading arm 
debris strea m at an observed peak distance of ^ 50 kpc (see 
Figure 6 of LMlOl) . However, such a rapidly moving satellite 
model yields radial velocities along the tidal debris streams 
that are systematically discrepant from observations by ^ 75 
km s"' . Similarly, it is neither possible to obtain a satisfac- 
tory fit for larger 8lsr values (©lsr ^ 280 km s"', for which 

The halo triaxiality is an added complication that has little bearing on 
the present discussion. 
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the halo mass scaling problem is even more extreme), nor for 
much lower values (0lsr 190 km s"', because the bary- 
onic bulge+disk mass component alone require 0lsr > 190 
km s"'). 

Another possibility we considered was to again fix the bary- 
onic mass (i.e., the bulge+disk component), but to a smaller 
value than the LMIO model, and allow the halo mass to vary. 
In particular, we attempted to fit a model with 8lsr = 232 km 
s"', but with the bulge+disk mass decreased by 10% from the 
ILMIO values. Even this small change in the baryonic mass 
required scaling up the dark matter halo mass by ^ 50% to 
keep 0LSR = 232 km s"' . The best-fit A^-body model in this 
case fit the Sgr trailing-tail velocities, but was a poor fit to the 
leading arm SDSS distances because of the increased speed of 
the Sgr core necessitated by the much larger halo. This illus- 
tration highlights the large changes in the halo mass effected 
by even small changes in the baryonic mass or the LSR ve- 
locity when fitting to observational data on the Sgr system. In 
fact, it is a testament to how well-constrained the Sgr system 
is by the current observational data that we are unable to fit 
the data if we change the dark halo model substantially. 

One possible way to construct an A^-body model of the Sgr 
dwarf that fits the observational data relatively well while dra- 
matically changing 9lsr is to adopt a Galactic halo whose 
scalelength is a factor of ^ 10 shorter than commonly adopted 
(from Thaio =12 kpc in the LMIO model to rhaio = 1 kpc). 
However, the Galactic rotation curve implied by such a short 
halo scalelength declines steeply outside t he solar circle (Fig - 
ure |20]|. in conflict with observations (e.g.. lSofue et al.ll2009l) . 
We therefore conclude that it is not possible to satisfacto- 
rily model the Sgr dwarf in a Milky Way model with the 
bulge and disk masses fixed at the LMIO values of Mbuige = 
3.4x lO'°M0,Mdisk= l.Ox lO"M0,andthe Milky Way halo 
scaled to produce Glsr much higher than 220 km s"^ Thus 
our (and other recent) suggestions that 8lsr is due an upward 
revision implies that the disk and/or bulge components - but 
not the halo - of the Milky Way are more massive than previ- 
ously thought. 

5. ABUNDANCES 

While the spectra in the Selected Areas were obtained pri- 
marily with kinematics in mind, they have sufficient resolu- 
tion and, for a large fraction of stars, sufficient S/N, to obtain 
information not only on metallicity but also abundance pat- 
terns. This allows us an independent estimate of abundance 
distributions th at, whil e of lo wer p recision than the echelle 
work of M onaco etal.1 (120051 120071) and iChou eTaP (l2007[ 
ITOlO). is derived for many Sgr stars, and is less biased than 
those M-giant studies. 

5.1. Sgr Metallicity 

Metallicities were measured for all stars using a software 
pipeline entitled "EZ_SPAM" (Easy Stellar Parameters and 
Metallicities), details of which will appear in a forthcom- 
ing paper (Carlin et al. 2011, in prep.). EZ_SPAM relies 
on the we ll-understood and ca l ibrated Lic k spectral indices 
(see, e.g.. IWorthey et al.l [19941 lFriellll987h to measure stel- 
lar properties from low-resolution spectra. In particular, esti- 
mates of [Fe/H] are derived for target stars using eight Fe in- 
dices combined with the H/3 index. Calibration of these multi- 
dimensional data comes from fits of known [Fe/H] values as 
a function of the Lick Fe and H/3 indices for stars in the atlas 
of .Schiavon (.2007> based on the spectra of ,Jones..l998.) . The 



EZ_SPAM code yields [Fe/H] measurements with la preci- 
sion of ^ 0.3 dex at S/N 20, decreasing to 0.1 dex at 
higher signal-to-noise (S/N > 50). 



Table 6 

Mean [Fe/H] for Sgr Debris in Kapteyn Selected Ai'eas of Tliis Study 
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<[Fe/H]> 


""[Fe/H] 


Npe/H]' 


Aq 
(degrees) 


71 


-1.14±0.19 


0.97±0.14 


24 


128.2 


94 


-1.13±0.08 


0.61±0.06 


57 


116.3 


93 


-1.25±0.1I 


0.47±0.08 


23 


103.2 


92 


-0.97±0.16 
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90.1 


117 


-1.08±0.08 


0.52±0.06 


43 


87.6 


116 


-1.21±0.2I 


0.64±0.15 


10 


74.9 



"Number of spectra with S/N > 20 providing reliably-measured [Fe/H]. 
''Cannot be measured for this field - too few spectra. 
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Figure 21. Measured values of [Fe/H] for all stars having spectra with S/N > 
20 in each SA field. The (red) hashed histogram is made up of Sgr candidates 
from our final samples in each field, and the solid black line represents all 
other stars (i.e., mostly Milky Way field stars, with perhaps some unidentified 
Sgr debris included) for which we obtained spectra. The distribution of Sgr 
metallicities is clearly different from that of the field stars in all of these 
regions (except perhaps SA 1 17, which is somewhat ambiguous), peaking at 
a more metal-poor mean value in each field. The bottom panel shows the 
metallicity distribution function (MDF) for all four fields in the 88° < Aq < 
128° portion of the trailing tail in our study. This fractional MDF consists of 
all the Sgr debris metallicities (red histograms) from the previous four panels, 
normalized by the total number of stars (147) in the sample. 

The metallicity distribution for all well-measured stars (i.e., 
those with spectra having S/N > 20) in each of the four fields 
(SAs 71, 94, 93, and 117) in which Sgr debris are reUably 
identified is given in Figure |2T| For each field, the solid 
(black) histogram shows [Fe/H] of non-Sgr stars, and the 
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hashed (red) histogram gives the distribution of [Fe/H] for 
stars selected to be Sgr members. The bottom panel repre- 
sents the distribution of metallicities for all Sgr members from 
the four trailing-tail fields, normalized by the total number of 
stars (147) in the sample to produce a fractional distribution. 
In each field, Sgr members are typically more metal-poor than 
the field stars, with the possible exception of those in SA 117. 

For each SA field in the survey, a maximum likelihood es- 
timate for [Fe/H] was derived from all well-measured stars in 
the final Sgr candidate sample. The resulting values for Sgr 
debris metallicities in each field are given in Table |6] along 
with cr[Fe/H]^ the dispersion in [Fe/H] about the mean. As 
was the case for the kinematics in SAs 92 and 116, we re- 
gard the [Fe/H] results in these fields (and, to a lesser de- 
gree, those in SA 71) with some skepticism, because the 
identification of Sgr debris in these fields is rather unreli- 
able. The mean metallicities for Sgr stars are displayed in 
Figure |22] as a function of A0; solid squares depict SAs 71, 
94, 93, and 117 (i.e., the "well-measured" fields), with open 
symbols included for SAs 92 and 116. Error bars repre- 
sent the uncertainties in the mean value from the maximum 
likelihood estimator; however, the scatter of [Fe/H] for Sgr 
candidates in each field is rather large. Typical fields have 
<j[Pe/H] = 0.5-0.6 dex about the quoted mean values, similar 
to the broad metallicity distribution function for Sgr sta rs seen 
by, e.g . , [Sm ecker-Hane & Mc William (2002), Zaggi a et al.l 
( 12004 . Isiegel et al., (,2007,) . and .Monaco et al., (.20051) . The 
scatter is even larger in SA 71 (at A0 = 128°); this may arise 
for a number of reasons. As can be discerned from Figure [1] 
SA 71 may be sampling Sgr debris stripped on multiple peri- 
centric passages (i.e., both gold and magenta debris may be 
present in this field). Furthermore, this is the lowest-latitude 
field among those in this study, and thus may be also suffering 
more contamination from Galactic thick disk stars. Finally, 
we note that S A 7 1 has been shown by Casetti-Dinescu et al. 
( 1200811) to contain a significant number of stars from the 
"Monoceros stream" overdensity, which could contribute to 
the inflation of the metallicity dispersion in this field, though 
it is unlikely that many Monoceros stars would lie within our 
Sgr radial velocity criteria for this field. Also shown in Fig- 
ure|22]is a solid line at constant [Fe/H] = -1.15, which is the 
mean value from the four well-measured fields; the tight cor- 
respondence of the mean values of each field to this line is 
consistent with the notion that debris along this narrow stretch 
of the trailing tail has constant metallicity. However, there is 
a hint of a shallow gradient, which we confirm by fitting a 
linear trend to the four good data points. This fit, overlaid 
as a dashed line in Figure |22] is [Fe/H] = -0.991±0.003 - 
(0.0014±0.0036) Aq. While suggestive of a slight gradient, 
the slope given is also consistent with zero within the errors 
of the fit. This is not surprising considering that nearly all de- 
bris in the portion of the stream contained within this study is 
expected to have been stripped on the same pericentric pas- 
sage of the Sgr core, as evidenced by the fact that all of our 
fields overlap gold-colored debris in Figure[T](i.e., debris that 
bec ame unbound during the last two perigalactic passages; 
see iLaw & Maiewskill2010al for more detail about the color 
scheme used). 

Our measured metallicity of [Fe/H] ~ -1.2 for Sgr trailing 
de bris is ^ 0.6 dex m ore metal-poor than the result obtained 
by lKelleretal.1 (12010b for traiUng-tail M-giants. This is not 
surprising, as M giants are biased toward relatively younger, 
more enriched stellar populations. In spite of this difference 
between the older, metal-poor Sgr stars in our SA fields and 



Fe/H] = -1.15 

■pe/H] = -0.9906 - 0.0014 Aq 




-1.5 - 



140 120 100 80 60 

Aq (degrees) 



Figure 22. Measured values of [Fe/H] for Sgr candidates in each SA field as 
a function of Sgr longitude, Aq . Filled squares and diamonds (and associ- 
ated error bars) show the maximum likelihood estimate from the individual 
Sgr candidates in each field (diamonds are the two fields lacking secure iden- 
tification of candidates). The mean value of the four well-measured fields 
(the filled squares), [Fe/H] = -1.15, is represented by the solid line that repro- 
duces the measurements well. A linear fit to the same four fields is shown as 
a dashed line, and is suggestive of a slight metallicity gradient of 1.4 X 10"' 
dex degree"' along the stream (though the fit is consistent with zero slope 
within the uncertainties). 

the more metal-rich M-giants, we measure a shallow gradient 
in [Fe/H] as a functio n of Aq, with a slope consistent with 
the lKeller et al.l (120 101) measurement, and just a simple offset 
in the zero-point metallicity. A more apt comparison for the 
m ean metallic i ty of Sgr debris in our SA fields is the work 
of iSesar et al] ([2010^, who used SDSS Stripe 82 data to de- 
velop a new technique for estimating metallicity from photo- 
metric data, which relies on combined information from both 
RR Lyrae variables and main-sequence stars from the same 
structure. Because SAs 94, 93, and 92 are within Stripe 82 
(and, in fact, we have used those SDSS data in our analysis), 
the Sesar et al. study probes identical stellar populations from 
the Sgr trailing tail as our work. This is borne out by the fact 
that our measured [Fe/H] = -1.15 is in excellent agreement 
wifli th e value of [Fe/H] = -1.20±0.1 derived by lSesaretall 
(12010') for Sgr debris along the same portion of the trailing 
stream. 

5.2. Metallicity Distribution Function 

Previous attempts to measure the metallicity distribution 
function (MDF) of the Sgr stream have suffered from both 
small number statistics and stellar t racers that h ave an inher- 
ent metallicity bias. For example, iChou et al.l (^007) mea- 
sured the MDF for M giant stars at several points along the 
Sgr leading stream as well as its core; using these data they 
attempted to reconstruct the MDF that the Sgr stream progen- 
itor would have had several Gyr ago. However, because M 
giants only form in metal rich stellar populations, this anal- 
ysis, while able to show that the mean metallicity of stars 
varies along the stream, was inadequate to assess the MDF 
across the full metallicity range of the system. In addition, 
the analysis was based on a relatively small sample (about 
7 dozen stars), nearly half of which are in the core of the 
Sgr dSph. Monaco et al. (2007) derived abundances of Sgr 
M giants along the trailing tidal tail, from which they derived 
([Fe/H]) = -0.61 ±0.13 between 80° < A© < 100° (from 6 
M giants), and ([Fe/H]) = -0.83 ± 0. 1 1 (mean of 4 stars) for 
debris even further along the trailing tail. The existence of 
a metallicity gradient among Sgr stream M giants along both 
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Figure 23. The normalized MDF from all 147 Sgr candidates in SAs 71, 94, 
93, and 1 17 with spectra having S/N > 20 is shown as a solid black histogram 
with grey fi ll. For comparison, in the upper panel we show the MDFs from 
IChou et aT l (2007) of the Sgr core (black dashed line) and Sgr leading arm 
M-giants (green dot-dashed Hne). In the lower panel the red (dot-dashed) 
and blue (dashed) histograms show the approximate Sgr MDF from several 
Gyr ago reconstructed by Chou et al. from linear combinations of the core 
and leading arm samples. The first of these was created by interpolating the 
MDFs at different orbital longitudes, and MDF 2 was created by assigning 
observed MDFs to particles in the , Law et al. 1 2005) Sgr model by the times 
they became unbound. 



the tra iling and leading tails was confirmed by Kell er et"an 
(120101) . who combined their additional measurements of 5 
stars at Aq = 66° (([Fe/H]) - -0.5) and 6 stars at A© = 132° 
(([Fe/H]) -0.7) with the Chou et al. and Monaco et al. re- 
sults to confirm the gradient in [Fe/H] among Sgr stream M 
giants. However, all of these M-giant studies suffer an inher- 
ent bias toward metal-rich stellar populations, and are likely 
not showing the true MDF of the Sgr system. 

Blue horizontal branch stars (BHBs) are another easily- 
identified and rather unambiguous tracer of halo substructure 
that has been used to probe the Sgr stream. However, BHB 
stars arise only in old, metal-poor populations, and are thus 
not ideal tracers of the global MDF of a system c onsisting of 
multiple stellar populations. Yannv et all (l2009bl) performed 
an extensive study of the Sgr tails using SDSS and SEGUE 
spectroscopy of BHB stars in both the northern and southern 
Galactic caps. BHB stars in the portions of both the leading 
(200° < Ao < 300°) and trailing (70° < Aq < 110°) tails in 
this study have MDFs peaking at ([Fe/H]} ^ -1.7, with sig- 
nificant numbers of stars as low as [Fe/H] ^ -2.5. Although 
this result turns up metal-poor populations not seen in M gi- 
ants, it is difficult to make conclusions about the overall MDF 
of the Sgr stream or progenitor based on biased metallicity 
tracers such as BHB stars and M giants. 

Fortunately, because the present analysis makes use of 
MSTO stars, it is far less susceptible to metallicity biases and 
can provide new insights into the MDF (particularly at the in- 
termediate to metal-poor end) of the stream (and therefore the 
progenitor) MDF. Of course, our spectra have ~ lOx worse 
resolution than the various echelle resolution studies, but our 
sample of Sgr stream stars is significantly larger, including 
147 with good enough S/N (> 20) for [Fe/H] measurements 
to the approximately < 0.3 dex level. As shown in Figure 1211 
and|22l the mean [Fe/H] for that portion of the stream probed 



by our S A data is about -1.1, but with a significant tail to both 
solar metallicity as well as very metal poor (< -2.0) stars. 
Indeed, our sample includes some rather metal-poor stars as- 
sociated with the Sgr system, with stars as metal-poor as the 
-2.5 dex BHBs seen by Y anny et all (12009b). 

Figure |23] shows the MDF derived from our data in the 
88° < Aq < 128° portion of the Sgr trailing tail. The MDF 
we derive is significantly broader and extending to much more 
metal-poor stars than indicated by the biased, M giant studies 
(shown in the upper panel of Figure |23] as a black dashed his- 
togram for the Sgr core and green dot-dashed lines for the 
leading arm), encompassing both the M- giant and BHB re- 
sults. A comparison of our MDF to the lChouetal] dlOOTh 
reconstruction of the Sgr M-giant MDF from several Gyr ago 
based on their core and leading-arm samples is given in the 
lower panel of Figure |23] Clearly our Sgr trailing-tail sam- 
ple is lacking the metal-rich component seen in the present- 
day core, but shows a similar distribution to the metal-poor 
tail of the reconstructed MDF. Additional metal-poor stars are 
present in our sample that are not seen in the M giant samples; 
these are likely draw n from similar populations to those in the 
lYann v et alj (I2009bl) study. 

Obviously, as has been shown by Chou et al. and others, the 
total MDF of the entire Sgr system will include a higher con- 
tribution of metal-rich stars when the core is included, but we 
also expect more metal-poor stars from those parts of the tails 
with larger separation from the core than we explore. Thus, 
while we cannot yet accurately reconstruct the total MDF of 
the Sgr system, at least we now have a better feel of the 
breadth of the Sgr MDF from the data shown in Figures |2T| 
and |23] Comparison of the latt er MDF with tho se of other 
MW dSphs (summarized, e.g., in lKirbv et al.l201 11) shows Sgr 
to be more typical of other MW satellites. In particular, the 
Sgr MDF resembles even more that of the LMC (as has been 
previo usl y suggested by, e.g., Monaco et al. 2003, Cole et d] 
I2005L and lMonaco et al. 2005), which has b een argued to be 
a chemical analog to the Sgr pr ogenitor b y Chou et al] ( 120101) 
and a morphological analog bv lLokas eta l. (2010). 

5.3. "Alpha" Abundances 

As shown in Section 5.1, we have identified metal-poor 
populations in (at least) four of the fields from our study, 
which explore a different segment of the stellar populations 
in the Sgr stream than previous M-giant studies. We have 
observed stars in these fields only at low resolution, and 
thus cannot do detailed element-by-element chemical analysis 
such as that enabled by high-resolution spectroscopy. How- 
ever, we can use the low-resolution Lick indices to explore 
relative a-abundances for the stars in our study. Specifically, 
we explore the relative Mg abundances using the Lick Mg b 
index centered at 5 160-5 190 A. Calibrating the Mg b index to 
an actual [Mg/Fe] abundance is difficult, because the strength 
of Mg lines is highly sensitive to surface gravity, with some 
additional sensitivity to effective temperature and [Fe/H]. Dis- 
entangling these effects is difficult with low-resolution spec- 
tra, but we can still explore a subset of the stars in our sam- 
ples in a way that is relatively free of the effects of surface 
gravity and temperature of individual stars. To do so, we se- 
lect only blue (0.2 <g-r < 0.7, or B-V < 0.9) stars, which 
should be mostly main-sequence dwarfs (thus, with similar 
surface gravity), since no giants are found at such blue col- 
ors. Furthermore, the temperature sensitivity of the Mg line 
strength, which is already much smaller than the log g sensi- 
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Figure 24. Relative values of Lick index ratio log (Mg b/<Fe>), where the indices are as described in the text and Carlin et al. 2011 (in prep ), for all blue 
stars (0.2 < ^j-r < 0.7 for SAs 94 and 93, and B-V < 0.9 for SAs 71 and 117) having spectra with5/A' > 30 in the four SA fields with securely identified Sgr 
debris. When comparing only predominantly dwarf stars of similar, blue photometric colors, the log (Mg b/<Fe>) ratio can be thought of as a proxy for [Mg/Fe], 
because variations in [Mg/Fe] with log g and color (i.e., temperature) are then minimized. Colored points represent all stars within the initial Sgr candidate RV 
selections, with color codes as in the legend. Black dots are all other stars outside the Sgr velocity range. For [Fe/H] > -1.5, Sgr candidates (colored points) 
typically occupy a region of lower Mg abundance at a given [Fe/H] than the black dots that are likely Galactic foreground stars. This behavior is typical for stars 
from most Galactic dSphs (relative to Galactic disk populations). At lower metallicities, the distributions converge. 



tivity, is mitigated by concentrating on a limited color range. 
In Figure|24]we show a "pseudo- [Mg/Fe]" ratio, given as the 
logarithm of the ratio of the Lick Mg b index to the mean of 
all eight Lick Fe indices (after transforming them to a com- 
mon scale), for all of the blue stars in SAs 71, 94, 93, and 
117 for which we have high enough signal-to-noise (> 30) to 
precisely measure indices and [Fe/H]. Black points in this di- 
agram are all stars with non-Sgr radial velocities, while our 
samples of all stars with Sgr-like RVs in each field are shown 
as colored points. It is readily apparent that the black (MW) 
points mostly occupy a different region of the diagram than 
the colored (Sgr) dots, which suggests an intrinsic chemical 
difference between the populations (though some overlap is 
expected, especially at low metallicities, where many MW 
halo stars likely resemble dSphs in their abundance patterns). 
Indeed, the behavior seen in Figure |24] is exactly that seen for 
many MW dSphs - for more metal-rich dSph stars, the Mg (or 
a) abundance is lower (on average) at a given [Fe/H] than in 
the Galactic populations, with the two populations converging 
at lower metallicities (i.e., at the "knee" in the dSph's distri- 
bution). Among the more metal-rich (and younger) M-giant 
populations of the Sgr stream, there is some indicati on that 
the knee in [g/Fe] oc curs at -1.2 < [Fe/H] < -1.0 ( Cho u et all 
I2010t [Monaco et al.i i2007). but this is difficult to assess be- 
cause of the lack of M-giants at lower metallicity. Thus the 
apparent convergence of Sgr trailing tail [Mg/Fe] with the 
plateau seen in Galactic stars at [Fe/H] < -1.5 may be an ex- 
tension of the same behavior seen in the M-giant studies. Al- 
ternatively, since we've already shown that the mean [Fe/H] 
along the trail ing tail differs between the M-giant sample of 
[Monaco et all (1^07), who find [Fe/H] -0.6, and our result 
of [Fe/H] ^ -1.2 (w hich is also consistent with the findings 
of ISesar et ai]l20Toh . our study may be sampling a distinctly 
older, more metal-poor population of Sgr debris than the M- 
giant tracers. Further characterization of the a-element be- 
havior along the Sgr trailing tail would benefit from either a 
calibration of our Mg b/<Fe> index onto [Mg/Fe] abundance 
or the identification of bona fide stream giant stars that are 
bright enough for echelle-resolution spectroscopic follow-up. 



6. SUMMARY 

We have presented the first large-scale study of the 3-D 
kinematics of the Sagittarius trailing tidal stream, with data 
spanning ~ 60° along the trailing tail. The data include deep, 
precise proper motions derived from photographic plates with 
a ~90-year baseline, and radial velocities from more than 
1500 low-resolution stellar spectra, of which > 150 have been 
identified as Sgr debris stars. Mean absolute proper motions 
of these Sgr stars in four of the six 40' x 40' fields from our 
survey have been derived with ^0.25-0.7 mas yr"' per field 
precision in each dimension (depending on the quality and 
depth of plate material and the number of spectra obtained 
in each field). Mean three-dimensional kinematics in each of 
these four fields have been shown to agree with the predicted 
Vgsr and /i;, from the Sagittarius disruption models of LMlOi 
However, there is a systematic disagreement in the fii cos(b) 
proper motions (with the exception of the somewhat problem- 
atical SA 71 field), which we use to assess refinements to the 
mass scale of the Milky Way (particularly its disk and bulge 
components). 

While proper motions along the portion of the trailing tail 
in this study provide constraints on Sgr tidal disruption mod- 
els, the fortuitous orientation of the Sgr plane also allows us 
to use the measured proper motions to derive the circular ve- 
locity at the Solar circle (or "Local Standard of Rest"), 8lsr- 
Our first-order approximation using only the /i/ cos(b) proper 
motions as constraints yields Glsr = 264 ±23 km s"'. From 
our measured 3-D kinematics, we find this fundamental Milky 
Way parameter to be 0lsr = 232 ± 14 km s~', or '--^ Icr higher 
than the lAU standard value of 220 km s~' . When we remove 
SA 71, a field in which it is more difficult to unambiguously 
identify Sgr debris, from the sample we find Glsr = 244 ±17 
km s"'. We suggest that the true value of Glsr lies some- 
where between 232-264 km s"' , while noting that all three of 
these estimates are consistent with each other within their Icr 
uncertainties. 

Our general result that the circular velocity at the Solar ra- 
dius is higher than the lAU standard of 220 km s"' agrees with 
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the recent derivation of 8lsr = 254 ±16 km s"' bv lReid et alj 
(|2P09) using trigonometric parallaxes of star forming regions 
in the outer disk. The sa me maser data from the Reid et al. 
study were reanalyzed by iBovy et a l. ( 2009), and yield a re- 
sult of 246 ± 30 km s"' (244 ± 13 km s"' if priors on the 
proper motion of Sgr A* are included, and 236 ±11 km s~' 
if the additional contribution of orbital fitting to the GD-1 
stellar stream is included). Again, our independent result is 
consistent with these studies, and inconsistent with the lAU 
accepted value of 220 km s~' for this fundamental constant 
at the 1-2(7 level. Identification of additional Sgr candidates 
could increase the accuracy of our determination of 0lsr, as 
could the addition of another epoch of accurate data to the 
proper motion measu rements. 

We note that while Reid et alj (l2009h argued that their mea- 
surement of 254 km s"' would imply an upward revision of 
the total Milky Way mass by a factor of ~ 2 (to a mass sim- 
ilar to that of M31), we have shown that this is not required 
to produce 0lsr even higher than that of Reid et al. Scal- 
ing the Milky Way dark matter halo up in mass to a level 
that yields 9lsr = 264 km s~' while simultaneously reproduc- 
ing known leading arm debris requires the Sgr core to have a 
high (^ 400 km s"') space velocity, resulting in RVs in the 
tidal streams that are discrepant by ^ 75 km s"' from mea- 
sured values. Instead, we show that because > 80% of the 
centripetal force at the location of the Sun is contributed by 
mass in the Galactic disk and bulge, an increase of ~ 50% 
in the mass of the disk-nbulge accounts for the additional ac- 
celeration needed to produce 264 km s~' rotation at the solar 
circle, while contributing only a small 7%) increase to the 
total virial mass of the Milky Way. With the additional con- 
straint on the disk-nbulge mass provided by our measurement 
of &LSR, we have found a satisfactory model of Sgr disruption 
that matches all of the constraints used in fitting the LMIO 
model, while additionally predicting a proper motion for the 
Sgr dwa rf that is in much better agreement with observations 
than the ILM Id model. 

Stellar metallicities have been derived from the low- 
resolution spectra of Sgr candidates, and the mean metallic- 
ity of Sgr tidal debris derived in each field. We find that a 
constant [Fe/H] = -1.15 is consistent with the observations of 
all four fields for which Sgr members were reliably identi- 
fied. However, a linear fit to these four data points suggests 
that a gradient of (1.4 x lO"-') dex degree"' is also reason- 
able (though this value is consistent with zero slope wi thin the 
uncertainty), in lin e with previous findings (e.g., Cho u et all 
120071: iKeller et al."2010) of a metallicity gradient among M- 
giants along both the leading and trailing tidal tails. The scat- 
ter of [Fe/H] in each of the survey fields is >0.5 dex, which 
is typical of the stellar populations seen in the cor e of the Sgr 
dSph (e.g., Smecker-Hane & McWilliam 2002; Zaggia et al.l 
120041; [Siegel et al.i i2007e iMonaco et al.. 2005). We show the 
metallicity distribution function for the trailing tail that is free 
from the biases inherent in previous studies of the Sgr MDF 
We find that the MDF of trailing debris is similar to MDFs 
of typical classical Milky Way dwarf spheroidals. A "pseudo- 
[Mg/Fe]" was measured based on the ratio of Lick Mg b and 
<Fe> indices; the behavior of log (Mg b/<Fe>) with [Fe/H] 
for Sgr main-sequence candidates is markedly different from 
Galactic stars of similar photometric colors (identified by ra- 
dial velocity) from among the same datasets. Furthermore, 
the trend is similar to that typically seen for dSphs, in that 
[Mg/Fe] is deficient at a given [Fe/H] for Sgr stars relative to 



the Milky Way field populations, and converges to a "knee" 
in Figure |24] at lower ([Fe/H] ^ -1.5) metallicity. This is a 
lower [Fe/H] than previously reported for Sgr M giants, and 
may reflect a bias intrinsic to those earlier M-giant studies. 
High-resolution spectroscopic follow-up will be necessary to 
confirm this trend among the old, metal-poor populations of 
recently-stripped Sgr debris. 

We appreciate the useful comments provided by the anony- 
mous referee. We thank Mei- Yin Chou for kindly sharing the 
Sgr MDF data used in Figure 23, and Heidi Newberg for many 
useful discussions. JLC acknowledges support from National 
Science Foundation grant AST-0937'523, and observing travel 
support from the NOAO thesis student program for proposal 
ID 2008B-0448. JLC and SRM acknowledge partial funding 
of this work from NSF grant AST-0807945 and NASA/JPL 
contract 1228235. DIC and TG acknowledge NSF grant AST- 
0406884. DRL acknowledges support provided by NASA 
through Hubble Fellowship grant # HF-5 1244.01 awarded by 
the Space Telescope Science Institute, which is operated by 
the Association of Universities for Research in Astronomy, 
Inc., for NASA, under contract NAS 5-26555. 

SRM is grateful to Allan Sandage for alerting him to the 
existence of the Mt. Wilson Selected Areas plates, for sug- 
gesting their use as a first epoch to expand our proper motion 
work in Kapteyn Selected Areas, and for his contribution of 
KPNO plates to our survey. 

Facilities: WIYN (Hydra), MMT (Hectospec), Sloan, 
Swope, MtW:1.5m, Du Pont, Mayall 

Funding for the SDSS and SDSS-II has been provided by 
the Alfred P. Sloan Foundation, the Participating Institutions, 
the National Science Foundation, the U.S. Department of En- 
ergy, the National Aeronautics and Space Administration, the 
Japanese Monbukagakusho, the Max Planck Society, and the 
Higher Education Funding Council for England. The SDSS 
Web Site is http://www.sdss.org/ 

The SDSS is managed by the Astrophysical Research Con- 
sortium for the Participating Institutions. The Participating 
Institutions are the American Museum of Natural History, As- 
trophysical Institute Potsdam, University of Basel, Univer- 
sity of Cambridge, Case Western Reserve University, Uni- 
versity of Chicago, Drexel University, Fermilab, the Institute 
for Advanced Study, the Japan Participation Group, Johns 
Hopkins University, the Joint Institute for Nuclear Astro- 
physics, the Kavli Institute for Particle Astrophysics and Cos- 
mology, the Korean Scientist Group, the Chinese Academy 
of Sciences (LAMOST), Los Alamos National Laboratory, 
flie Max-Planck-Institute for Asti'onomy (MPIA), the Max- 
Planck-Institute for Astrophysics (MPA), New Mexico State 
University, Ohio State University, University of Pittsburgh, 
University of Portsmouth, Princeton University, the United 
States Naval Observatory, and the University of Washington. 

REFERENCES 

Abadi, M. G., NavaiTO, J. P., Steinmetz, M., & Eke, V. R. 2003, ApJ, 591, 
499 

Abazajian, K. N., Adelman-McCarthy, J. K., Agiieros, M. A., et al. 2009, 

ApJS, 182, 543 
Alard, C. 2001, A&A, 377, 389 
Andrae, R. 2010, ai'Xiv: 1009.27551 

Bedin, L. R., Piotto, G., King, 1. R., & Anderson, J. 2003, AJ, 126, 247 

Bellazzini, M., Ferraro, R R., & Ibata, R. 2003, AJ, 125, 188 

Bellazzini, M., Newberg, H. J., Corrend, M., Ferraro, F. R.. & Monaco, L. 

2006, A&A, 457, L21 
Belokurov, V., Zucker, D. B., Evans, N. W., et al. 2006, ApJ, 642, L137 
Belokurov, V., Evans, N. W., twin, M. J., et al. 2007, ApJ, 658, 337 



26 



Carlin et al. 



Bovy, J., Hogg, D. W., & Rix, H. 2009, ApJ, 704, 1704 

Bullock, J. S., & Johnston, K. V. 2005, ApJ, 635, 931 

Carlin, J. L., Casetti-Dinescu, D. L, Grillmair, C. J., Majewski, S. R., & 

Girard, T. M. 2010, ApJ, 725, 2290 
Carraro, G., Zinn, R., & Moni Bidin, C. 2007, A&A, 466, 181 
Casetti-Dinescu, D. I., Carlin, J. L., Girard, T. M., et al. 2008, AJ, 135, 2013 
Casetti-Dinescu, D. I., Girard, T. M., Majewski, S. R., et al. 2009, ApJ, 701, 

L29 

Casetti-Dinescu, D. I., Majewski, S. R., Girard, T. M., et al. 2006, AJ, 132, 
2082 

Chou, M., Cunha, K., Majewski, S. R., et al. 2010, ApJ, 708, 1290 
Chou, M., Majewski, S. R., Cunha, K., et al. 2007, ApJ, 670, 346 
Cole, A. A., Tolstoy, E., Gallagher, 111, J. S., & Smecker-Hane, T. A. 2005, 
AJ, 129, 1465 

Correnti, M., Bellazzini, M., Ibata, R. A., Ferraro, F. R., & Varghese, A. 

2010, ApJ, 721, 329 
Dehnen, W., & Binney, J. J. 1998, MNRAS, 298, 387 
Dinescu, D. I., Girard, T. M., van Altena, W. F, & Lopez, C. E. 2005, ApJ, 

618, L25 

Dinescu, D. 1., Majewski, S. R., Girard, T. M., & Cudworth, K. M. 2000, AJ, 
120, 1892 

Dinescu, D. I., Majewski, S. R., Girard, T. M., et al. 2002, ApJ, 575, L67 

Dohm-Palmer, R. C, Helmi, A. Morrison, H. L., et al. 2001, ApJ, 555, L37 

Fabricant, D., Fata, R., Roll, J., et al. 2005, PASP, 1 17, 141 1 

Feast, M., & Whitelock, P 1997, MNRAS, 291, 683 

Fellhauer, M., Belokurov, V., Evans, N. W., et al. 2006, ApJ, 651, 167 

Font, A. S., Johnston, K. V., Bullock, J. S., & Robertson, B. E. 2006, ApJ, 

638, 585 
Friel, E. D. 1987, AJ, 93, 1388 

Frinchaboy, P. M., Muiioz, R. R., Phelps, R. L., Majewski, S. R., & Kunkel, 

W. E. 2006, AJ, 131,922 
Geisler, D., Smith, V. V., Wallerstein, G., Gonzalez, G., & Charbonnel, C. 

2005, AJ, 129, 1428 
Ghez, A. M., Salim, S., Weinberg, N. N., et al. 2008, ApJ, 689, 1044 
Girardi, L., Grebel, E. K., Odenkirchen, M., & Chiosi, C. 2004, A&A, 422, 

205 

Grillmair, C. J. 2006a, ApJ, 645, L37 
— . 2006b, ApJ, 651, L29 
— . 2009, ApJ, 693, 1118 

Grillmair, C. J., & Dionatos, O. 2006, ApJ, 643, L17 

Hargreaves, J. C, Gilmore, G., Irwin, M. J., & Carter, D. 1994, MNRAS, 

269, 957 
Helmi, A. 2004, ApJ, 610, L97 

Helmi, A., & White, S. D. M. 2001, MNRAS, 323, 529 
Ibata, R., Lewis, G. F, Irwin, M., Totten, E., & Quinn, T. 2001, ApJ, 551, 
294 

Ibata, R. A., Gilmore, G., & Irwin, M. J. 1994, Nature, 370, 194 

Ibata, R. A., Irwin, M. J., Lewis, G. F, Ferguson, A. M. N., & Tanvir, N. 

2003, MNRAS, 340, L21 
Johnston, K. V., Law, D. R., & Majewski, S. R. 2005, ApJ, 619, 800 
Johnston, K. V., Spergel, D. N., & Hernquist, L. 1995, ApJ, 451, 598 
Johnston, K. V., Zhao, H., Spergel, D. N., & Hernquist, L. 1999, ApJ, 512, 

L109 

Jones, L. A. 1998, PhD thesis. The University of North Carolina at Chapel 
Hill 

Kalirai, J. S., Richer, H. B., Hansen, B. M., et al. 2004, ApJ, 601, 277 
Kapteyn, J. C. 1906, Plan of selected areas (Groningen, Hoitsema brothers, 
1906.) 

Keller, S. C, Yong, D., & Da Costa, G. S. 2010, ApJ, 720, 940 

Kerr, F J., & Lynden-Bell, D. 1986, MNRAS, 221, 1023 

Kirby, E. N., Lanfranchi, G. A., Simon, J. D., Cohen, J. G., & Guhathakurta, 

P 2011, ApJ, 727,78 
Kleyna, J., Wilkinson, M. 1., Evans, N. W., Gilmore, G., & Frayn, C. 2002, 

MNRAS, 330, 792 
Koposov, S. E., Rix, H., & Hogg, D. W. 2010, ApJ, 712, 260 
Kuijken, K., & Tremaine, S. 1994, ApJ, 421, 178 
Law, D. R., Johnston, K. V., & Majewski, S. R. 2005, ApJ, 619, 807 
Law, D. R., & Majewski, S. R. 2010a, ApJ, 714, 229 (LMIO) 
— . 2010b, ApJ, 718, 1128 

Law, D. R., Majewski, S. R., & Johnston, K. V. 2009, ApJ, 703, L67 
Lokas, E. L., Kazantzidis, S., Majewski, S. R., et al. 2010, ApJ, 725, 1516 



Majewski, S. R. 1992, ApJS, 78, 87 
— . 1993, ARA&A, 31, 575 

Majewski, S. R. 1999, in Astronomical Society of the Pacific Conference 
Series, Vol. 165, The Third Stromlo Symposium: The Galactic Halo, ed. 
B. K. Gibson, R. S. Axelrod, & M. E. Putman, 76-H 

Majewski, S. R., Law, D. R., Polak, A. A., & Patterson, R. J. 2006, ApJ, 
637, L25 (MLPP) 

Majewski, S. R., Munn, J. A., & Hawley, S. L. 1996, ApJ, 459, L73-I- 

Majewski, S. R., Skrutskie, M. F, Weinberg, M. D., & Ostheimer, J. C. 

2003, ApJ, 599, 1082 

Majewski, S. R., Kunkel, W. E., Law, D. R., et al. 2004, AJ, 128, 245 
Martmez-Delgado, D., Gomez-Flechoso, M. A., Apaiicio, A., & Carrera, R. 

2004, ApJ, 601,242 

Mihalas, D., & Binney, J. 1981, Galactic astronomy: Structure and 

kinematics /2nd edition/ (W. H. Freeman and Co.) 
Mink, D. J., Wyatt, W. F, Caldwell, N., Comoy, M. A., Furesz, G., & 

Tokarz, S. P. 2007, in Astronomical Society of the Pacific Conference 

Series, Vol. 376, Astronomical Data Analysis Software and Systems XVI, 

ed. R. A. Shaw, F Hill, & D. J. Bell, 249-+ 
Miyamoto, M., & Nagai, R. 1975, PASJ, 27, 533 
Monaco, L., Bellazzini, M., Bonifacio, R, et al. 2007, A&A, 464, 201 
Monaco, L., Bellazzini, M., Bonifacio, R, et al. 2005, A&A, 441, 141 
Monaco, L., Bellazzini, M., Ferraro, F. R., & Pancino, E. 2003, ApJ, 597, 

L25 

Muiioz, R. R., Carlin, J. L., Frinchaboy, P M., et al. 2006a, ApJ, 650, L51 
Muiioz, R. R., Majewski, S. R., & Johnston, K. V. 2008, ApJ, 679, 346 
Muiioz, R. R., Majewski, S. R., Zaggia, S., et al. 2006b, ApJ, 649, 201 
Munn, J. A., Monet, D. G., Levine, S. E., et al. 2004, AJ, 127, 3034 
— . 2008, AJ, 136, 895 

Newberg, H. J., Yanny, B., Rockosi, C, et al. 2002, ApJ, 569, 245 
Oiling, R. P, & Merrifield, M. R. 1998, MNRAS, 297, 943 
Peiiarrubia, J., Belokurov, V, Evans, N. W., et al. 2010, MNRAS, 408, L26 
Pryor, C, & Meylan, G. 1993, in Astronomical Society of the Pacific 

Conference Series, Vol. 50, Structure and Dynamics of Globular Clusters, 

ed. S. G. Djorgovski & G. Meylan, 357 
Pryor, C, Piatek, S., & Olszewski, E. W. 2010, AJ, 139, 839 
Reid, M. J., & Brunthaler, A. 2004, ApJ, 616, 872 
Reid, M. J., Menten, K. M., Zheng, X. W., et al. 2009, ApJ, 700, 137 
Robin, A. C, Reyle, C, Derriere, S., & Picaud, S. 2003, A&A, 409, 523 
Sbordone, L., Bonifacio, P, Buonanno, R., et al. 2007, A&A, 465, 815 
Schiavon, R. R 2007, ApJS, 171, 146 

Schlegel, D. J., Finkbeiner, D. P, & Davis, M. 1998, ApJ, 500, 525 
Scares, F. H., Kapteyn, J. C, van Rhijn, P. J., Joyner, M. C, & Richmond, 
M. L. 1930, Mount Wilson catalogue of photographic magnitudes in 
selected areas 1-139 (Carnegie institution of Washington) 
Searle, L., & Zinn, R. 1978, ApJ, 225, 357 
Sesar, B., Ivezic, Z. Grammer, S. H., et al. 2010, ApJ, 708, 717 
Shetrone, M., Venn, K. A., Tolstoy, E., et al. 2003, AJ, 125, 684 
Shetrone, M. D., Cote, P, & Sargent, W. L. W. 2001, ApJ, 548, 592 
Siegel, M. H., Dotter, A., Majewski, S. R., et al. 2007, ApJ, 667, L57 
Siegel, M. H., et al. 2011, ApJ, accepted (arXiv:l 108.6276 1 
Smecker-Hane, T. A., & McWilliam, A. 2002, ArXiv:020541 1 
Sofue, Y, Honma, M., & Omodaka, T. 2009, PASJ, 61, 227 
Sohn, S. T., Majewski, S. R., Muiioz, R. R., et al. 2007, ApJ, 663, 960 
Tolstoy, E., Hill, V., & Tosi, M. 2009, ARA&A, 47, 371 
Tolstoy, E., Venn, K. A., Shetrone, M., et al. 2003, AJ, 125, 707 
Tonry, J., & Davis, M. 1979, AJ, 84, 1511 
van Leeuwen, F. 2007, A&A, 474, 653 

Venn, K. A., Irwin, M., Shetrone, M. D., et al. 2004, AJ, 128, 1177 
Vogt, S. S., Mateo, M., Olszewski, E. W., & Keane, M. J. 1995, AJ, 109, 151 
Worthey, G., Faber, S. M., Gonzalez, J. J., & Burstein, D. 1994, ApJS, 94, 
687 

Yanny, B., Newberg, H. J., Grebel, E. K., et al. 2003, ApJ, 588, 824 
Yanny, B., Rockosi, C, Newberg, H. J., et al. 2009a, AJ, 137, 4377 
Yanny, B., Newberg, H. J., Johnson, J. A., et al. 2009b, ApJ, 700, 1282 
Yuan, F, Zhu, Z., & Kong, D. 2008, Chinese J. Astron. Astrophys., 8, 714 
Zaggia, S., Bonifacio, P., Bellazzini, M., et al. 2004, Memorie della Societa 
Astronomica Italiana Supplementi, 5, 291 



